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1.  INTRODUCTION 

During  the  past  decade*  the  study  of  artificially  layered  structures  has  become  an  increasingly 
important  subject  in  condensed  matter  physics  and  solid  state  electronics.  Progress  in  the 
controlled  growth  of  semiconductor  heterostructures  and  compositionally  modulated  metals 
has  resulted  in  the  synthesis  of  many  new  types  of  physical  structures.  Experimental  and 
theoretical  studies  of  inversion  layers  at  semiconductor  interfaces  have  provided  a  better 
understanding  of  the  physical  properties  of  confined  two-dimensional  electron  gases.  These 
studies  have  also  led  to  the  discovery  of  new  and  interesting  phenomena  such  as  the 
fractional  quantum  Hall  effect,  and  to  the  invention  and  practical  utilization  of  novel 
electronic  devices  such  as  the  high  electron  mobility  GaAs  transistor.  Future  progress  in 
VLSI  and  submicron  electronics  will  undoubtedly  reflect  the  results  of  current  research  on 
semiconductor  heterostructures,  passivating  oxide  overlayers,  metal-semiconductor 
interfaces,  and  metallic  multilayers. 

For  the  past  several  years  we  have  been  studying  the  electronic  structure  and  related  physical 
properties  of  interfaces  that  play  a  key  role  in  microelectronics.  We  have  been  particularly 
interested  in  understanding  the  relationship  between  the  atomic-scale  structure  of  interfaces 
and  the  local  electronic  structure,  such  as  the  energy  level  spectrum  and  the  spatial  range  of 
interface  states.  Since  we  are  dealing  with  very  complex  systems,  it  is  necessary  to 
represent  actual  interfaces  by  highly  simplified  but  instructive  atomic-scale  models.  Using 
these  models,  we  can  carry  out  first-principles  electronic  structure  calculations,  obtaining 
detailed  information  about  localized  interface  states  and  related  optical  and  magnetic 
properties.  In  many  cases  our  theoretical  results  can  be  used  directly  to  understand  and 
interpret  specific  experimental  results,  or  indirectly,  to  predict  general  trends. 

During  the  period  covered  by  this  contract,  we  studied  interfaces  between  (a)  different 
semiconductors,  such  as  Ge  and  GaAs;  (b)  semiconductors  and  insulators,  such  as  Si  and 
Si02;  (c)  metals  and  semiconductors,  such  as  Pd2Si  and  Si;  and  (d)  different  metals,  such  as 
Co  and  Cr.  We  also  examined  interfaces  between  crystalline  and  amorphous  tetrahedrally 
coordinated  semiconductors.  Our  earlier  studies  were  concerned  primarily  with  electronic 
and  optical  properties,  while  more  recent  ones  focused  on  electronic  and  magnetic  properties. 

In  the  next  several  sections  we  will  briefly  summarize  our  major  research  studies,  including 
the  objectives,  the  general  approach,  the  principal  conclusions,  and  where  the  subject  stands 
today.  References  to  our  publications  are  given  at  the  end  of  the  first  section  in  which  they 
are  mentioned. 

Additional  information  may  be  found  in  the  representative  research  papers  and  reviews  that 
are  reproduced  in  the  Appendices.  We  have  not  included  all  the  papers  we  have  published 
during  this  period,  but  only  those  that  we  feel  will  be  of  interest  to  the  readers  of  this  Final 
Report.  We  have  excluded  highly  theoretical  papers  as  well  as  those  already  published  in 
widely  available  scientific  journals  such  as  the  Journal  of  Chemical  Physics,  Physical  Review, 
or  Journal  of  Vacuum  Science  and  Technology.  Some  of  the  papers  appearing  in  the 
Appendices  have  been  published  in  conference  proceedings  and  research  journals  that  are  not 
as  widely  available,  while  others  are  still  in  press. 
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2.  ELECTRONIC  STRUCTURE  OF  {110!  Ge/GaAs  INTERFACES 

Our  earliest  study  in  this  field1  was  devoted  to  the  prototypic  atomically  abrupt 
he tero junction  between  lattice-matched  semiconductors,  namely,  the  non-polar  (110) 
interface  between  Ge  and  GaAs.  For  simplicity,  we  considered  [U0|  Ge/GaAs  superlattices 
having  repeat  periods  of  8  and  12  layers.  There  are  two  atoms  per  layer,  so  that  the  unit 
cells  contain  1$  and  24  atoms.  Of  these,  8  atoms  lie  at  the  surfaces  of  the  Ge  and  GaAs 
slabs,  while  the  remaining  (8  and  Id)  atoms  lie  inside  these  slabs.  The  8-layer  model  is  the 
simplest  that  can  be  used  to  study  [110]  Ge/GaAs  interfaces,  while  the  12-layer  model  is 
somewhat  more  accurate.  Both  models  are  pictured  in  Fig.  1  below. 

If  one  calculates  the  electronic  structure  of  these  superlattices,  one  would  expect  to  find 
localized  interface  states  associated  with  the  Ge-Ga  and  Ge-As  bonds.  Such  Interface  states 
do  indeed  exist,  but  we  find  that  none  lies  inside  the  range  of  the  thermal  gap,  so  that  none 
of  these  states  is  electronically  active,  acting  as  a  trap  or  as  a  recombination  center.  This  is 
our  most  significant  result.  Moreover,  we  find  that  all  atomic  layers  are  electrically  neutral 
(there  is  no  interlayer  charge  transfer),  and  that  the  interface  orbitals  are  quite  compact, 
being  confined  almost  entirely  to  the.  interfacial  layers.  For  this  reason  adjacent  interfaces 
in  our  superlattice  models  are  effectively  isolated  from  one  another.  In  short,  the 
discontinuities  in  chemical  composition  (Ge  to  GaAs)  and  crystal  structure  (diamond  to 
sphalerite)  are  not  sufficiently  large  to  produce  Interface  states  in  the  thermal  gap. 


16  Atoms/Unit  Cell 


Fig.  1  Atomic  arrangements  of  [110]  Ge/GaAs  superlattices  with  16  and  24  atoms  per  unit 
cell  The  numbers  to  the  right  and  below  the  atoms  indicate  the  coordinates  in  the  direction 
perpendicular  to  the  page.  In  the  lower  panel,  the  bulk  Ge  layers  are  not  shown  but  their 
arrangement  is  readily  inferred.  After.  Ref.  1. 
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Further  studies  suggest  that  interface  states  in  the  thermal  gap  can  only  be  produced  by 
substantially  larger  discontinuities,  such  as  those  associated  with  interfacial  vacancies  or 
misfit  dislocations  and  characterized  by  dangling  bonds. 

Since  the  band  gaps  of  Ge  and  GaAs  differ  by  several  tenths  of  an  eV,  It  would  be  useful  to 
know  how  much  of  this  band  offset  takes  place  at  the  conduction  band  edge,  and  how  much 
at  the  valence  band  edge.  Numerical  results  for  band  offsets  are  very  sensitive  to  the 
computational  details  of  the  calculation.  Our  calculations1  were  not  sufficiently  detailed  to 
allow  us  to  obtain  accurate  values  for  the  band  offsets,  but  qualitatively  we  found  that  the 
conduction  band  offset  is  considerably  larger  than  the  valence  band  offset.  So  far  as  we  are 
aware,  accurate  numerical  results  for  band  offsets  have  not  yet  been  obtained  for  any 
heterostructure  by  first-principles  methods,  nor  have  systematically  accurate  results  been 
obtained  for  a  large  class  of  heterostructures  by  semi-empirical  methods. 

Although  our  own  work  and  that  of  others  has  elucidated  the  nature  of  Idealized  (1101 
Ge/GaAs  interfaces,  the  experimental  situation  has  proved  to  be  far  more  complicated  than 
anybody  imagined  originally.  Because  of  interatomic  diffusion  and  interfacial  chemical 
reactions,  the  atomic-scale  structure  of  actual  Ge/GaAs  interfaces  is  still  not  understood. 
Band  offsets  have  been  measured  for  various  crystallographic  orientations,  bnt  the  measured 
results  are  undoubtedly  influenced  by  structural  and  chemical  abnormalities  at  the  interfaces 
that  are  ignored  in  the  idealized  models  used  by  ourselves  and  others.  It  is  difficult  to  make 
further  theoretical  progress  on  Ge/GaAs  heterostructures  because  of  the  uncertain  nature 
of  the  atomic-scale  structure. 

There  have  been  many  experimental  attempts  to  grow  atomically  abrupt  Ge/GaAs  Interfaces, 
but  interfacial  diffusion  and  interfacial  chemical  reactions  have  frustrated  all  of  these 
efforts.  The  Ge/GaAs  system  continues  to  be  of  great  scientific  interest  because  of  its 
inherent  simplicity  —  at  least  in  principle  —  but  this  system  has  not  yet  been  exploited  for 
practical  purposes  because  of  the  difficulties  encountered  in  growing  atomically  abrupt, 
structurally  perfect  interfaces.  Most  practical  research  focuses  on  the  more  complicated 
GaAs/GaxAlj.xAs  system,  which  is  more  easily  grown  In  near-perfect  form. 

Further  information  as  well  as  detailed  references  concerning  Ge/GaAs  heterostructures 
may  be  found  in  our  original  paper,1  in  our  summary  of  PCSI-9*,  and  in  two  recent  review 
articles.3'4  Note  that  Ref.  4  Is  included  in  this  report  as  Appendix  A. 

REFERENCES 

1.  F.  Herman  and  R.V.  Kasowski,  "Electronic  Structure  of  (1101  Ge-GaAs  Superlattices 

and  Hetero junctions,"  Phys.  Rev.  B17,  472  (1978). 

2.  F.  Herman,  "Recent  Developments  in  the  Physics  and  Chemistry  of  Semiconductor 

Interfaces,"  in  Proceeding  PCSt-9  (American  Vacuum  Society,  New  York,  1982),  p. 
xli. 

3.  F.  Herman,  "Physics  and  Chemistry  of  Semiconductor  Interfaces:  Some  New 

Directions,"  J.  Vac.  ScL  Techn.  21,  643  (1982). 

4.  F.  Herman,  "Electronic  Structure  of  Semiconductor  Interfaces,"  J.  Phys.  (Paris) 

Colloque  45  C5-375  (1984)  (cf.  Appendix  A). 
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3.  ELECTRONIC  STRUCTURE  OF  Si/SI02  INTERFACES 

The  interface  between  Si  and  its  native  oxide  is  undoubtedly  the  most  important  interface  in 
modern  solid  state  electronics.  Passivating  oxide  layers  protect  integrated  circuits  from 
environmental  contaminants  and  also  serve  to  isolate  individual  circuits  from  neighboring 
ones.  Even  though  there  is  a  vast  literature  on  the  Si/Si02  interface,  our  understanding  of 
this  interface  is  still  rudimentary.  Although  the  density  of  localized  interface  states  has  been 
reduced  empirically  to  an  acceptable  minimum  for  present-day  applications,  there  is  still 
considerable  interest  in  understanding  the  structural  and  chemical  features  responsible  for 
residual  interface  states.  It  is  anticipated  that  such  knowledge  will  he  used  eventually  to 
control  or  totally  eliminate  residual  interface  states  in  future  device  applications,  for  which 
present  levels  would  be  intolerable. 

In  order  to  come  to  grips  with  this  interface,  we  developed  an  Idealized  structural  model 
which  ignores  structural  disorder  In  the  oxide  and  replaces  the  amorphous  SIO,  with  a 
crystalline  form  of  Si02  (cristobalite)  that  can  be  lattice-matched  to  crystalline  Si.*“®  In  this 
way  we  establish  a  computationally  tractable  model  for  the  interface  between  crystalline  Si 
and  crystalline  Si02.  According  to  this  model,  roughly  half  the  interfacial  Si  atoms  are 
bonded  to  the  Si02  region,  while  the  remaining  half  are  not.  Because  of  the  open  nature  of 
SI02,  It  is  postulated  that  H,  O,  and  OH  groups  can  move  through  the  oxide,  and,  upon 
reaching  the  Interface,  saturate  most  of  the  dangling  bonds  arising  from  the  partially  bonded 
Si  atoms  at  the  interface.  Clearly,  the  motion  and  concentration  of  these  species  will  depend 
strongly  on  growth  conditions. 

On  the  basis  of  this  simple  model  of  the  Si/Si02  interface,  we  investigated  the  nature  of 
localized  electronic  states  arising  from  various  types  of  structural  defects  located  at  or  near 
this  interface.6-8  The  results  are  too  detailed  to  summarize  here,  but  some  defects  were 
shown  to  give  rise  to  levels  that  lie  within  the  SI  thermal  gap,  while  others  produced  levels 
that  overlap  the  Si  valence  and  conduction  bands  and  hence  are  electrically  inactive.  The 
most  important  electronically  active  defects  are  unsaturated  Si  atoms  at  the  interface.  In 
order  to  make  further  progress,  it  would  be  necessary  to  deal  not  only  with  the  Si  substrate 
and  the  oxide  overlayer,  but  also  with  the  H,  O,  and  OH  species  that  take  part  in  the 
formation  of  the  interface. 

Recent  experimental  work  has  shown  that  the  SI/S102  interface  is  almost  atomically  abrupt, 
spanning  one  or  two  atomic  layers  of  strained  Si  substrate  and  one  or  two  Si02  layers.  In 
our  idealized  model,  the  interface  is  atomically  abrupt,  but  there  is  considerable  lattice  strain 
at  the  interface  which,  if  relieved,  would  extend  the  effective  range  of  the  interface  to  a  few 
atomic  layers  on  either  side,  in  agreement  with  experimental  indications. 

In  more  recent  work  (cf.  Section  6),  we  made  a  start  toward  constructing  interfaces  between 
crystalline  and  amorphous  semiconductors.  Judging  by  the  results  obtained  thus  far  for  the 
interface  between  crystalline  and  amorphous  Si,  we  expect  to  be  able  to  obtain  an 
atomic-scale  model  for  the  Interface  between  crystalline  Si  and  amorphous  Si02.  This  model 
would  be  more  realistic  than  the  idealized  model  described  above,  in  that  structural  disorder 
at  the  interface  Is  taken  Into  account  explicitly.  However,  we  would  stili  expect  a  significant 
number  of  dangling  Si  bonds  to  be  present  at  or  near  the  interface.  We  believe  that  these 
dangling  bonds  are  the  principal  source  of  localized  interface  states  in  the  thermal  gap  of  Si. 
In  practice,  we  would  expect  more  of  these  states  to  be  eliminated  by  the  saturation  of 
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dangling  bonds  by  H,  O,  and  OH  groups  than  by  structural  rearrangements  (lattice 
relaxation).  We  believe  that  further  insight  into  the  nature  of  the  Si/Si02  interface  will 
come  only  after  the  formation  and  growth  of  this  interface  is  modeied  theoretically,  taking 
into  account  the  role  played  by  co-existing  H,  O,  and  OH  species. 

For  further  information  and  detailed  references  concerning  the  Si/Si02  interface,  see  Refs. 
5  through  8.  Note  that  Ref.  7  appears  in  this  report  as  Appendix  B. 
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(cf.  Appendix  B). 

8.  F.  Herman  and  R.V.  Kasowski,  "Electronic  Structure  of  Defects  at  Si-Si02  Interfaces," 

J.  Vac.  ScL  Techn.  19,  395  (1981). 


4.  IMPROVED  COMPUTATIONAL  METHODS  FOR  COMPLEX  CRYSTALS 

In  addition  to  preparing  a  review  paper  on  computational  methods  for  surfaces,  overlayers, 
and  interfaces,0  we  wrote  two  original  theoretical  papers  suggesting  improvements  in  the 
treatment  of  complex  crystals  by  the  extended  muffin-tin  orbital  (EMTO)  method.10,11  By 
complex  crystals  we  mean  not  only  bulk  crystals  containing  a  large  number  of  atoms  per  unit 
cell,  but  also  superlattices,  including  atomic  layers  chemisorbed  on  crystalline  substrates  and 
represented  by  repeating  slab  geometries.  At  present,  self-consistent  EMTO  programs  are 
being  refined  (principally  by  Kasowski  and  his  associates)  so  that  these  programs  could  be 
used  to  determine  total  energies  and  optimal  atomic  geometries.  We  expect  the 
improvements  proposed  in  Refs.  10  and  11  to  be  incorporated  in  next-generation  EMTO 
computer  programs. 

The  detailed  formalism  of  our  theoretical  approach  appears  in  Ref.  10.  The  essential  ideas 
are  summarized  in  Ref.  1 1,  which  is  included  in  this  report  as  Appendix  C. 

9.  F.  Herman,  "Electronic  Structure  Calculations  of  Interfaces  and  Overlayers  in  the 

1980’s,"  J.  Vac.  Sci.  Techn.  16,  1101  (1979). 

10.  F.  Casula  and  F.  Herman,  "Generalized  Muffin-Tin  Orbitals  for  Electronic  Structure 

Studies  of  Surfaces,  Interfaces,  and  Organic  Solids,"  J.  Chem.  Phys.  78,  858  (1983). 

11.  F.  Herman  and  F.  Casula,  "Electronic  Structure  Calculations  for  Complex  Crystals," 

Nuovo  Cimento  2D,  1636  (1983)  (cf.  Appendix  C). 
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5.  ELECTRONIC  STRUCTURE  AND  SCHOTTKY  BARRIER  HEIGHTS  AT 
Pd2Si/Si(lll)  INTERFACES 

One  of  the  most  interesting  problems  in  the  physics  of  semiconductor  interfaces  is  the 
nature  of  rectifying  metal-semiconductor  contacts.  In  particular,  what  are  the  mechanisms 
primarily  responsible  for  rectifying  (Schottky)  barriers  at  such  contacts?  In  order  to  come 
to  grips  with  this  subject,  we  decided  to  study  a  particularly  simple  but  technologically 
important  system,  namely,  epitaxially  grown  Pd  overlayers  on  Si(lll)  substrates.4*12 

On  the  basis  of  LMTO  electronic  structure  calculations,  we  concluded  that  hybridized  Pd-Si 
states  at  the  interface  could  not  by  themselves  account  for  measured  Schottky  barriers 
heights,  and  that  another  mechanism  or  structural  feature  had  to  be  invoked.  After 
examining  various  possibilities,  we  conjectured  that  Pd  atoms  moved  into  the  Si  substrate  in 
advance  of  a  growing  interface,  and  that  these  Pd  atoms  could  form  triangular  clusters,  these 
being  precursors  to  the  arrangement  of  Pd  atoms  in  Pd2Sl.  According  to  our  model,  two 
adjacent  Si  atoms  are  replaced  by  a  triangular  Pd  cluster  lying  in  the  plane  bisecting  the  bond 
between  the  two  Si  atoms  that  have  been  removed.  That  is  to  say.  Si  divacancies  are  replaced 
by  Pd3  clusters. 

We  also  noted  that  Pd3  clusters  are  more  easily  accommodated  in  hexagonal  Si  than  in  cubic 
Si  because  of  the  different  shapes  of  the  interstitial  cavities  in  these  two  structures.  We 
reasoned  that  the  presence  of  Pd3  clusters  in  the  region  ahead  of  the  growing  interface 
would  cause  the  Si  layers  to  grow  with  the  hexagonal  rather  than  the  cubic  stacking 
sequence.  The  net  effect  would  be  a  transition  layer  of  hexagonal  Si  between  the  Pd2Si 
overlayer  and  the  cubic  Si  substrate.  Since  hexagonal  Si  has  a  smaller  band  gap  than  cubic 
Si,  the  effective  band  gap  at  the  Pd2/hexagonal-Si  interface  would  be  a  few  tenths  of  an  eV 
less  than  at  a  Pd2/cubic-Si  interface,  this  reduction  being  sufficient  to  account  for  the 
measured  Schottky  barrier  height  for  this  system. 

Some  attempts  have  been  made  to  observe  triangular  Pd  clusters  and  hexagonal  stacking 
experimentally,  but  so  far  there  is  no  direct  evidence  in  favor  of  this  model.  However,  the 
idea  of  stacking  disorder  has  been  used  recently  to  account  for  some  features  of  the 
reconstructed  Si(lU)  7x7  surface.  Interstitial  Pd  atoms  have  also  been  invoked  in  other 
models  of  the  Pd2Si/Si(lll)  interface. 

The  epitaxially  grown  Pd2Si/Si(ll  1)  interface  continues  to  be  of  interest  because  of  the 
nearly  perfect  lattice  match,  and  because  only  one  silicide  forms  when  Pd  is  reacted  with  Si. 
Many  different  models  for  the  interfacial  atomic  structure  have  been  proposed,  but  there  is 
still  no  common  agreement  among  different  investigators.  It  is  anticipated  that  much  more 
work  will  have  to  be  done  before  this  system  — »  which  is  one  of  the  simplest  —  is 
understood.  We  believe  that  our  own  work  is  significant  because  it  is  the  first  to  call 
attention  to  the  possibility  of  stacking  disorder  in  the  neighborhood  of  internal  surfaces  as  a 
possible  factor  determining  Schottky  barrier  heights. 

For  further  details,  see  our  original  paper,12  which  appears  in  this  report  as  Appendix  D. 
See  also  Ref.  4,  which  is  reproduced  in  Appendix  A. 

12.  F.  Herman,  F.  Casuia,  and  R.Y.  Kasowski,  "Electronic  States  and  Schottky  Barriers  at 
Pd2Sl/Sl(lll)  Interfaces,"  Physics  117B-118B,  837  (1983)  (cf.  Appendix  D). 
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6.  INTERFACES  BETWEEN  CRYSTALLINE  AND  AMORPHOUS  TETRAHEDRALLY 
COORDINATED  SEMICONDUCTORS 

Ball-and-stlck  models  of  disordered  regions  and  of  interfaces  between  ordered  and 
disordered  regions  hare  been  constructed  by  hand  by  various  investigators.  We  do  not  regard 
these  models  as  trustworthy  because  they  represent  the  outcome  of  a  long  series  of 
subjective  and  largely  arbitrary  human  judgments.  One  can  usually  fashion  a  structure  that 
satisfies  one’s  own  preconceived  notions.  For  example,  a  model  of  amorphous  Si  containing 
no  dangling  bonds  can  be  built  in  the  form  of  a  continuous  random  network,  but  is  this  what 
actually  occurs  in  nature?  We  thought  it  would  be  useful  to  develop  a  fully  automatic  method 
for  constructing  interfaces  between  crystalline  and  amorphous  semiconductors,  thereby 
avoiding  subjective  influences.  Our  approach  is  based  on  a  set  of  rules  according  to  which 
such  interfaces  can  be  constructed  by  computer  without  human  intervention.  Of  course,  *' 
rules  embodied  in  the  computer  programs  represent  plausible  physical  and  chev  J 
constraints,  such  as  tetrahedral  coordination,  most  likely  bond  lengths  and  bond  angles,  e 

For  mathematical  purposes  we  begin  with  a  superlattice  composed  of  alternating  slat 
ordered  and  disordered  close-packed  hard  spheres.  The  interaction  between  spheres  is 
represented  by  a  Lennard-Jones  potential.  The  ordered  slabs  have  the  face-centered  cubic 
structure,  while  the  disordered  slabs  are  characterized  by  random  dose-packing.  These  slabs 
may  be  regarded  as  co-existing  frozen  and  melted  regions  of  hard-sphere  systems  held 
together  by  Lennard-Jones  forces. 

Treating  the  centers  of  these  spheres  as  the  vertices  of  regular  or  irregular  tetrahedra,  we 
can  select  a  subset  of  these  centers  which  forms  a  network  of  tetrahedra  and  in  contact  only 
at  their  vertices.  That  is  to  say,  adjacent  tetrahedra  touch  at  vertices  but  not  along  common 
edges.  In  the  ordered  regions  this  tetrahedral  network  (or  simplidal  graph)  includes  all  the 
sphere  centers,  while  in  the  disordered  regions  most  but  not  all  the  sphere  centers  belong  to 
the  tetrahedral  network.  We  now  examine  each  of  the  tetrahedra  in  turn  and  discard  those 
having  highly  irregular  shapes.  We  then  place  Si  atoms  at  the  vertices  of  all  the  tetrahedra 
still  belonging  to  the  network.  We  also  place  atoms  at  the  centroids  of  these  tetrahedra. 
Finally,  we  join  each  centroid  atom  to  its  four  vertex  atoms,  thus  establishing  tetrahedrally 
coordinated  chemical  bonds  between  centroid  and  neighboring  vertex  atoms. 

In  the  ordered  regions,  this  procedure  generates  a  diamond  lattice,  with  each  centroid  atom 
and  each  vertex  atom  lying  at  the  center  of  a  regular  tetrahedron.  In  the  disordered  regions 
and  at  interfaces,  this  procedure  generates  a  discontinuous  random  network  which  is  for  the 
most  part  tetrahedrally  coordinated.  Here  all  centroid  atoms  have  four  nearest  neighbors 
lying  at  the  vertices  of  their  tetrahedra.  However,  not  all  vertex  atoms  ..re  shared  by  four 
touching  tetrahedra,  so  that  many  vertex  atoms  have  fewer  than  four  nearest  neighbor 
atoms.  At  the  completion  of  this  first  stage  of  the  construction,  the  atomic  positions  have 
not  yet  been  relaxed,  still  being  determined  by  the  frozen  and  melted  hard  sphere 
Lennard-Jones  svstems  with  which  we  started.  We  have  already  published  a  brief  account  of 
this  first  stage. 

At  the  end  of  Stage  I,  the  interface  extends  over  a  few  atomic  layers.  The  average 
coordination  number  in  the  interior  of  disordered  regions  is  3,  and  this  rapidly  increases  to  4 
in  the  neighborhood  of  interfaces.  Thus,  the  model  is  characterized  by  many  dangling  bonds 
in  the  disordered  regions.  The  number  of  dangling  bonds  decreases  rapidly  as  one  passes 


Research  Summary 


9 


through  the  interfacial  region  into  the  ordered  regions.  We  have  also  calculated  the 
electronic  density  of  states  for  this  model  using  the  recursion  method,  and  we  find  a  direct 
correlation  between  the  number  of  dangling  bonds  in  a  given  layer  and  the  layer  density  of 
electronic  states  in  the  thermal  gap. 

The  second  stage  of  the  construction  is  now  in  progress.  During  this  stage  the  model  is 
improved  in  two  important  respects.  First,  the  lattice  is  relaxed  in  accordance  with  a 
semi-empirical  potential  energy  function  of  the  Keating  type.  Second,  bonds  are  formed 
between  nearby  pairs  of  atoms  which  have  dangling  bonds  at  the  end  of  Stage  I.  Most  but 
not  all  of  the  Stage  I  dangling  bonds  are  eliminated  during  Stage  II,  leading  to  a  considerably 
more  realistic  atomic-scale  model.  In  the  third  stage,  now  being  designed,  there  will  be  a 
random  interchange  of  adjacent  vertex  and  centroid  atoms  in  the  disordered  regions  and  at 
the  interfaces,  so  that  the  original  distinction  between  these  two  types  of  atoms  will  be 
eliminated.  The  interchange  is  being  patterned  after  that  used  recently  by  Weaire  and 
Wooten  in  their  computer  generated  models  of  bulk  amorphous  Si. 

The  construction  as  it  stands  applies  equally  to  interfaces  between  crystalline  and  amorphous 
C,  Si,  Ge,  and  grey  Sn.  The  construction  is  readily  generalized  to  deal  with  interfaces 
between  crystalline  Si  and  amorphous  Si02.  The  model  automatically  includes  dangling  bonds 
at  and  near  the  interfaces.  One  can  saturate  some  or  all  of  these  dangling  bonds  by  diffusing 
H,  O,  and  OH  groups  into  the  structure  in  accordance  with  prescribed  statistical  rules.  This 
can  be  done  during  or  after  the  second  stage  of  construction,  thereby  simulating  different 
oxide  and  interface  growth  conditions. 

We  believe  that  computer  modeling  of  the  growth  of  semiconductor  interfaces  will  provide 
new  insights  into  the  nature  of  the  structural  defects  that  occur  naturally  during  interface 
growth.4  The  present  study,13  which  appears  here  as  Appendix  E,  represents  a  step  In  this 
direction. 

13.  F.  Herman  and  P.  Lambin,  "Electronic  Structure  of  Interfaces  between  Crystalline  and 
Amorphous  Tetrahedrally  Coordinated  Semiconductors,"  Proceedings  of  the  17th 
International  Conference  on  the  Physics  of  Semiconductors  (Springer- Verlag,  Berlin,  1984) 
in  press  (cf.  Appendix  E). 
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7.  FERROMAGNETIC-ANTIFERROMAGNETIC  INTERFACES  AND  MULTILAYERS 

During  the  past  year  we  broadened  our  interests  to  include  interfacial  problems  involving 
magnetic  materials.  As  a  means  for  getting  into  this  new  field  and  learning  the  experimental 
situation,  we  wrote  a  review  paper  dealing  with  exchange  and  magnetocrystalline 
anisotropy.14  During  a  search  for  idealized  models,  our  attention  was  drawn  to  jellium,  which 
can  be  paramagnetic,  ferromagnetic,  or  antiferromagnetic  depending  on  the  electron  density. 
This  led  us  to  study  phase  transitions  in  jellium.1®  Our  work  yielded  more  precise  estimates 
of  the  densities  at  which  phase  transitions  occur,  but  also  indicated  that  the  stability  range 
for  the  ferromagnetic  phase  could  be  very  narrow  or  nonexistent  altogether.  We  had 
originally  intended  to  construct  interfaces  between  ferromagnetic  and  antiferromagnetic 
jellium,  but  since  our  calculations  indicated  that  the  stability  range  of  the  ferromagnetic 
state  of  jellium  was  too  uncertain  for  our  purposes,  we  abandoned  this  idea. 

We  then  decided  to  study  simplified  but  realistic  interfaces  rather  than  highly  idealized 
electron  gas  models.  Since  we  wanted  to  study  dose-packed  metals,  we  switched  from  the 
extended  muffin-tin  orbital  method  —  which  is  particularly  useful  for  open  structures  —  to 
the  LMTO/ASA  method  of  O.K.  Andersen  —  which  is  among  the  most  powerful  methods 
presently  available  for  carrying  out  first-principles,  spin-polarized,  self-consistent  electronic 
and  magnetic  structure  calculations  for  complex  metallic  systems. 

Our  first  set  of  LMTO/ASA  calculations  dealt  with  interfaces  between  ferromagnetic  NIjFe 
and  antiferromagnetic  MnFe.1®  To  avoid  problems  associated  with  structural  disorder,  we 
represented  alloys  by  suitably  chosen  ordered  compounds,  and  on  this  basis  studied  the 
exchange  coupling  across  bilayers  composed  of  Permalloy  (75  percent  Ni  and  25  percent  Fe) 
and  50-50  MnFe  alloy.  In  practice,  we  studied  spin  distributions  in  superlattices  composed  of 
the  ordered  compounds  Ni3Fe  (Cu3Au  structure)  and  MnFe  (rocksalt  structure).  The  spin 
distributions  that  we  obtained  are  consistent  with  experiment. 

Our  second  set  of  LMTO/ASA  calculations  dealt  with  interfaces  between  ferromagnetic  Co 
and  antiferromagnetic  Cr.17,1**  We  were  again  interested  in  understanding  exchange 
coupling  and  spin  distributions  near  ferromagnetic/antiferromagnetic  interfaces.  The  Cr 
regions  were  represented  by  a  theoretically  stabilized  commensurate  antiferromagnetic  spin 
arrangement  which  simulates  its  experimental  counterpart,  the  impurity-stabilized 
commensurate  antiferromagnetic  phase  of  Cr.  In  the  case  of  single  Co  layers  separated  by 
Cr  multilayers,  the  coupling  across  the  interfaces  was  ferromagnetic,  while  for  Co 
multilayers  separated  by  Cr  multilayers,  the  coupling  could  be  ferromagnetic  or 
antiferromagnetic.  The  spin  arrangements  in  the  Cr  regions  could  be  weakly  or  strongly 
antiferromagnetic  depending  on  the  crystallographic  orientation  of  the  interfaces,  and  on  the 
degree  to  which  the  commensurate  antiferromagnetic  phase  of  Cr  is  stabilized  by  the 
addition  of  suitable  transition  metal  impurities. 

We  believe  that  Co/Cr  multilayers  will  eventually  be  synthesized  in  the  laboratory,  and  our 
calculations  are  intended  to  provide  some  insight  into  the  magnetic  structure  of  such 
multilayers,  as  a  guide  to  the  design  and  interpretation  of  future  experiments.  For  further 
information,  see  Refs.  17  and  18.  Ref.  17,  an  abbreviated  version  of  Ref  18,  appears  here 
as  Appendix  F.  We  are  also  in  the  process  of  writing  two  additional  papers  on  magnetic 
problems  associated  with  interfaces.19,20 
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14.  N.H.  March,  P.  Lambin,  and  F.  Herman,  "Exchange  and  Magnetocrystalline  Anisotropy 

in  3d  Cooperative  Magnets  in  Single*  and  Two-Phase  Metallic  Systems,"  J.  Magnetism 
and  Magnetic  Materials  (in  press) 

15.  F.  Herman  and  N.H.  March,  "Cooperative  Magnetism  in  Metallic  Jellium  and  in  the 

Insulating  Wigner  Electron  Crystal,"  Solid  State  Comm.  50,  725  (1984). 

16.  P.  Lambin  and  F.  Herman,  "Electronic  and  Magnetic  Structure  of  Idealized  Metallic 

Multilayers:  Ni3Fe/FeMn  System,”  Phys.  Rev.  B  (in  press). 

17.  F.  Herman,  P.  Lambin,  and  O.  Jepsen,  "Electronic  Structure  of  Idealized 
Cobalt-Chromium  Multilayers,"  J.  Appl.  Phys.  (in  press):  Proc.  30th  Annual  Conf. 
Magnetism  and  Magnetic  Materials  (cf.  Appendix  F). 

18.  F.  Herman,  P.  Lambin,  and  O.  Jepsen,  "Electronic  Structure  of  Idealized 
Cobalt-Chromium  Superlattices,"  submitted  to  Phys.  Rev.  B 

19.  F.  Herman,  P.  Lambin,  and  N.H.  March,  "Response  of  Antiferromagnetic  Metals  to 

Changes  in  Exchange  Field"  (in  preparation). 

20.  F.  Herman,  P.  Lambin,  and  N.H.  March,  "Description  of  Magnetic  Inhomogeneities  in 

Ferromagnetic  Metals"  (in  preparation). 


8.  OTHER  PUBLICATIONS 

Carburized  tungsten  is  of  considerable  interest  because  of  its  catalytic  properties.  In  order 
to  get  some  insight  into  the  nature  of  this  material,  experimental  efforts  were  undertaken  at 
Stanford  to  understand  the  electronic  structure  of  crystalline  tungsten  carbide,  which 
provides  a  convenient  structural  model  for  carburized  tungsten.  As  part  of  an  IBM-Stanford 
collaboration,  we  calculated  the  energy  band  structure  of  tungsten  carbide  and  contributed  to 
the  analysis  of  angular  resolved  photoemission  spectra  of  tungsten  carbide  single  crystals. 
Many  of  the  salient  features  of  these  spectra  were  successfully  interpreted  on  the  basis  of 
our  theoretical  band  structure.21*22  We  also  contributed  to  a  theoretical  study  of  the 
electronic  structure  of  idealized  lattice-matched  Ge/Al  interfaces.23  We  published  the 
outline  of  two  survey  lectures  on  semiconductor  interfaces  presented  at  the  1983  Brazilian 
Summer  School  on  Semiconductor  Physics,24  as  well  as  a  historical  account  of  our  early  days 
in  semiconductors23  based  on  an  interview  conducted  by  Dr.  Kris  Szymborski  on  behalf  of 
the  International  Project  on  the  History  of  Solid  State  Physics. 


21.  P.M.  Stefan,  M.L.  Shek,  I.  Lindau,  W.E.  Spicer,  F.  Herman,  L.I.  Johansson,  and  G. 

Brogren,  "Electronic  Structure  and  Photoemission  Studies  of  Tungsten  Carbide,"  J. 
Vac.  Sci.  Techn.  20,  622  (1982). 

22.  P.M.  Stefan,  M.L.  Shek,  I.  Lindau,  W.E.  Spicer,  L.I.  Johansson,  F.  Herman,  R.V. 

Kasowski,  and  G.  Brogren,  "A  Photoemission  Study  of  Tungsten  Carbide  (0001)," 
Phys.  Rev.  B  24,  5423  (1984). 

23.  I.P.  Batra  and  F.  Herman,  "Electronic  States  in  Model  Ge-Al  Interfaces,"  J.  Vac.  Sci. 

Techn.  Al,  1080  (1983). 

24.  F.  Herman,  "Electronic  Structure  of  Semiconductor  Interfaces,"  Brazilian  J.  Phys., 

Special  Volume:  Physics  of  Semiconductors,  (1983),  p.  335. 

25.  F.  Herman,  "Elephants  and  Mahouts  —  Early  Days  in  Semiconductor  Physics,"  Physics 

Today  37,  56,  (1984). 
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9.  PROFESSIONAL  ACTIVITIES 

Member,  Organizing  Committee,  Workshop  on  Effective  One-Electron  Potentials  in  Real 
Materials,  Ossining,  New  York,  Mar.  21-22,  1980 
Member,  Organizing  Committee,  Annual  Conferences  on  the  Physics  and  Chemistry  of 
Semiconductor  Interfaces  (PCSI),  1980-1984 

Conference  Chairman,  Ninth  Annual  Conference  on  the  Physics  and  Chemistry  of 
Semiconductor  Interfaces  (PCSI-9),  Asiiomar,  California,  Jan.  27-29,  1982 
Lecturer,  Summer  School  of  Semiconductor  Physics,  University  of  Campinas,  Brazil,  Feb. 
8-9,  1983 

Member,  Executive  Committee,  17th  International  Conference  on  the  Physics  of  Semi¬ 
conductors,  San  Francisco,  California,  Aug.  6-10,  1984 


10.  PRESENTATIONS 

Electronic  Structure  of  [110\  Ge-GaAs  Superlattices  and  Interfaces 

•  Sth  Annual  Conference  on  the  Physics  of  Compound  Semiconductor  Interfaces,  UCLA, 
Los  Angeles,  California,  Jan.  17,  1978 

•  American  Physical  Society,  Washington,  D.C.,  Mar.  29,  1978 
Electronic  Structure  of  Model  Si-SiOj  Interfaces 

•  International  Topical  Conference  on  the  Physics  of  Si02  and  its  Interfaces,  IBM 
Yorktown  Heights,  New  York,  Mar.  22,  1978  (with  I.P.  Batra  and  R.V.  Kasowski) 

•  Navy  Workshop  on  Surfaces  and  Interfaces,  Rockwell  Science  Center,  Thousand  Oak, 
California,  Nov.  15,  1979 

•  Workshop  on  Si-Si02  Interfaces,  Stanford  University,  Dec.  13,  1979 

•  American  Physical  Society,  New  York,  Mar.  26,  1980  (with  D.J.  Henderson  and  R.V. 
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Electronic  Structure  of  Cuprous  Chloride 

•  American  Physical  Society,  Chicago,  Mar.  23,  1979  (with  R.V.  Kasowski) 

Critique  of  Xa  and  Xafi  Methods 

•  Workshop  on  Effective  One-Electron  Potentials  in  Real  Materials,  Ossining,  New 
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Electronic  Structure  of  Vacancies  and  Interstitials  in  5/0 2 
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•  American  Physical  Society,  Phoenix,  Arizona,  Mar.  19,  1981 
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Interpretation  of  Angie-Integrated  and  Angle-Resolved  Photoemission  Measurements  on  Tungsten 
Carbide 

•  American  Physical  Society,  Phoenix,  Arizona,  Mar.  19,  1981  (with  P.  Stefan  and  W.E. 
Spicer) 

Electronic  Structure  of  the  Si-SiOj  Interface 

•  Condensed  Matter  Group,  Cambridge  University,  United  Kingdom,  Apr.  22,  1981 

•  Max-Planck-Institute  for  Solid  State  Research,  Stuttgart,  Germany,  Apr.  24,  1981 

•  International  Conference  on  Insulating  Films  on  Semiconductors,  Erlangen,  Germany, 
Apr.  27,  1981  (Invited  Paper) 

•  University  of  Dortmund,  Germany,  May  1,  1981 

•  Laboratory  for  the  Physics  of  Solids,  University  of  Paris  VI,  France,  May  5,  1981 

•  Institute  of  Solid  State  Physics,  Autonomous  University,  Cantobianco,  Madrid,  Spain, 
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Angle-Resolved  Photoemission  Results  from  Tungsten  Carbide  \0001\, 
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W.E.  Spicer) 

Electronic  Bam /  Structure  of  Tungsten  Carbide 

•  American  Physical  Society,  Dallas,  Texas,  Mar.  8,  1982  (with  R.V.  Kasowski,  P. 
Stefan,  and  W.E.  Spicer) 

Recent  Developments  in  the  Physics  and  Chemistry  of  Semiconductor  Interfaces 

•  Surface  Science  and  Catalysis  Science  Seminar,  Lawrence  Radiation  Laboratory, 
Berkeley,  California,  Mar.  18,  1982 

•  Solid  State  Physics  Seminar,  University  of  California  at  Irvine,  Apr.  5,  1982 

•  Solid  State  Physics  Seminar,  Institute  of  Physics,  Univ.  of  Cagliari,  Italy,  Sept.  20, 
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Electronic  Structure  and  Schottky  Barriers  at  Pd 2  Si-Si(lll)  Interfaces 

•  16th  International  Conference  on  the  Physics  of  Semiconductors,  Montpellier,  France, 
Sept.  8,  1982  (with  F.  Casula  and  R.V.  Kasowski) 

•  American  Physical  Society,  Los  Angeles,  California,  Mar.  21,  1983 

Electronic  Structure  Calculations  for  Complex  Crystals 

•  5th  International  Conference  on  Ternary  and  Multinary  Compounds,  Cagliari,  Italy, 
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•  Laboratory  for  Scientific  Computation,  Rio  de  Janeiro,  Brazil,  Feb.  2,  1983 

Two  Lectures  on  Semiconductor  Interfaces 

•  Summer  School  on  Semiconductor  Physics,  University  of  Campinas,  Brazil,  Feb.  8  and 
9,  1983 

Electronic  Structure  of  Semiconductor  Interfaces 

•  Internationa]  Conference  on  the  Dynamics  of  Interfaces,  Lille,  France,  Sept.  15,  1983 
(Invited  Paper) 

•  Solid  State  Physics  Seminar,  Max  Planck  Institute  for  Solid  State  Research,  Stuttgart, 
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Cooperative  Magnetism  in  Metallic  Jellium  and  in  the  Insulating  Wlgner  Electron  Crystal 

•  American  Physical  Society,  Detroit,  Mar.  26,  1984  (with  N.H.  March) 

Electronic  and  Magnetic  Structure  of  Idealised  Ferromagnetic  -  Antiferromagnetic  Multilayers 

•  American  Physical  Society,  Detroit,  Mar.  26,  1984  (with  P.  Lambin) 
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ELECTRONIC  STRUCTURE  OF  SEMICONDUCTOR  INTERFACES* 

Frank  Herman,  IBM  Research  Laboratory 
San  Jose,  California  95193  USA 


Rdsumd.  —  Bien  que  nous  ayons  aujourd’hui  une  bonne  comprehension  de  la  structure 
dlectronique  des  interfaces  iddaux  entre  semi-conducteurs,  incluant  la  nature  des  6 tats 
localises  d’interface,  des  etudes  experimentales  rdcentes  suggerent  que  la  structure  atomique 
des  interfaces  reels  est  considerablement  plus  complexe  que  celle  decrite  par  les  modiles 
idealises  utilises  en  recherche  theorique  actuelle.  Pour  comprendre  plus  complhtement  les 
proprietes  electroniques  des  interfaces,  11  est  done  essentiel  de  ddvelopper  une  modilisation 
plus  realists  de  la  structure  atomique  de  ces  interfaces.  Nous  illnstrerons  les  points  faibles  et 
les  points  forts  des  theories  d’aujourd'hui  par  une  discussion  de  la  structure  electronique  des 
interfaces  Ge/GaAs,  Si/Si02  et  Pd2Sl/Si(lll).  Nous  discuterons  egalement  le 
developpement  de  modiles  atomiques  amellorees  d’interface  entre  semi-conducteurs  cristal- 
lins  et  amorphes. 


Abstract  —  Although  we  now  have  a  good  understanding  of  the  electronic  structure  of 
idealized  semiconductor  interfaces,  including  the  nature  of  localized  interface  states,  recent 
experimental  studies  suggest  that  the  atomic  structure  of  actual  interfaces  is  considerably 
more  complicated  than  that  described  by  the  idealized  models  used  in  current  theoretical 
research.  If  we  are  to  understand  the  electronic  properties  of  Interfaces  more  fully,  it  is 
essential  that  we  develop  more  realistic  structural  models  for  interfaces.  We  will  illustrate 
the  strengths  and  weaknesses  of  present-day  theories  by  discussing  the  electronic  structure 
of  Ge/GaAs,  Si/Si02,  and  Pd2Si/Si(Ul)  interfaces.  We  will  also  discuss  the  development 
of  improved  atomic-scale  models  for  interfaces  between  crystalline  and  amorphous  semicon¬ 
ductors. 


I.  INTRODUCTION 


One  of  the  most  exciting  and  rapidly  expanding  subjects  in  solid  state  physics  is  the  study  of 
semiconductor  heterostructures  and  superlattices.  4  Improved  crystal  growth  techniques 
such  as  molecular  beam  epitaxy,8  and  the  development  of  refined  surface-sensitive  spectros¬ 
copies,  have  stimulated  wide-ranging  scientific1”8  and  technological4  research  efforts  on 


such  systems  as  Ge/GaAs,  AlxGaj_xAs/GaAs,  and  InAs/GaSb.  Recent  improvements  in 
sample  preparation  and  measurement  techniques  have  also  inspired  fundamental  studies  of 
technologically  important  interfaces  such  as  oxide-passivated  silicon  surfaces6'8  and 
rectifying  metal-semiconductor  contacts.9,10 


In  this  paper  we  will  sketch  the  present  state  of  our  understanding  of  the  electronic  structure 
of  some  important  types  of  semiconductor  interfaces.  Most  of  our  current  ideas  are  based 
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on  simple,  highly  idealized  atomic  interface  models.  For  purposes  of  discussion,  we  will 
define  the  ideal  semiconductor  interface  as  the  atomically  abrupt  boundary  between  two 
lattice-matched,  crystallographically  compatible  constituents.  It  is  assumed  that  each  of  the 
constituents  has  its  bulk  crystal  structure  right  up  to  the  interface,  that  there  is  exact  lattice 
registry  at  the  interface,  and  that  there  are  no  dangling  bonds  at  or  near  the  interface.  If 
the  two  constituents  are  not  exactly  lattice-matched,  it  is  assumed  that  the  residual  strain  is 
relieved  by  well-separated  misfit  dislocations.  By  using  crystals  or  alloys  having  different 
chemical  compositions,  crystal  structures,  and  doping  profiles,  it  is  possible  to  synthesize 
heterostructures  having  a  wide  range  of  physical  characteristics. 

Recent  experimental  studies  indicate  that  some  semiconductor  systems,  such  as 
AlxGaj.xAs/GaAs1  and  strained-layer  superlattices,11  can  be  made  smooth  on  the  scale  of  a 
single  atomic  layer,  and  can  be  grown  almost  free  of  imperfections.  The  availability  of 
nearly  ideal  interfaces  and  interfacial  arrays  (heterostructures  and  superlattices)  is  of  course 
of  the  greatest  scientific  and  technological  importance.  Most  interfaces,  however,  contain  a 
wide  variety  of  structural  and  chemical  imperfections;  much  current  research  is  concerned 
with  the  identification,  modification,  or  reduction  of  these  imperfections. 


II.  SEMICONDUCTOR  HETEROJUNCTIONS:  Ge/GaAs  INTERFACES 

Because  of  its  apparent  simplicity,  the  Ge/GaAs  system  has  been  widely  studied.  Ge  and 
GaAs  have  compatible  crystal  structures  and  nearly  identical  lattice  constants  (Ge:  5.65735 
A;  GaAs:  5.6537  A).  Since  the  band  gaps  of  GaAs  and  Ge  differ  by  0.9  eV,  and  the 
conduction  band  edges  occur  at  different  positions  in  the  reduced  zone,  there  must  be  a  band 
structure  discontinuity  at  the  interface  which  can  be  represented  by  valence  and  conduction 
band  offsets.  In  principle,  these  offsets  can  be  different  for  interfaces  having  different 
crystallographic  orientations. 

It  is  important  to  know  the  band  offsets  because  they  determine  the  depth  of  quantum  wells 
in  heterostructures  such  as  GaAs/Ge/GaAs  and  in  superlattices.  Theoretical  estimates12*13 
of  band  offsets  and  experimental  measurements14  are  only  in  rough  agreement  with  one 
another.  This  is  hardly  surprising:  some  of  the  estimates  are  based  on  oversimplified 
models,12  while  the  more  detailed  estimates13  are  usually  based  on  idealized  interface  models 
which  may  differ  in  some  important  respects  from  real  interfaces. 

The  discontinuity  in  crystal  potential  associated  with  the  changes  in  crystal  structure  and 
chemical  composition  across  the  interface  may  give  rise  to  localized  electronic  states.  When 
the  energies  of  such  states  lie  within  the  thermal  gap,  they  can  serve  as  trapping  and 
recombination  centers,  and  hence  are  of  considerable  importance.  (By  thermal  gap  we  mean 
the  energy  range  between  the  highest  occupied  valence  band  level  and  the  lowest  unoccupied 
conduction  band  level)  Since  the  Ge/GaAs  interface  discontinuity  is  substantially  weaker 
than  the  discontinuity  at  a  free  Ge  or  GaAs  surface,  interface  states  are  pushed  out  of  the 
valence  and  conduction  bands  to  a  much  lesser  extent  than  intrinsic  surface  states  are. 

Theoretical  studies  of  ideal  non-poiar  (110)  Ge/GaAs  Interfaces13*17  demonstrate  that 
localized  interface  states  do  not  extend  into  the  thermal  gap.  But  localized  states  associated 
with  the  Ge-Ga  and  Ge-As  bonds  at  the  Interface  do  occur  within  the  "stomachs"  of  the 
projected  valence  and  conduction  bonds  and  in  the  forbidden  band  just  outside  thermal  gap. 
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Localized  interface  states  associated  with  Ge  overlayers  on  GaAs  have  been  detected  by 
photoemission.18 

In  the  case  of  ideal  polar  (100)  Ge/GaAs  interfaces,  theoretical  studies18*20  indicate  that 
localized  interface  states  form  a  band  that  extends  into  the  thermal  gap.  Since  the  Fermi 
level  passes  through  this  interface  band,  this  ideal  interface  would  exhibit  metallic  behavior, 
contrary  to  experimental  results.  It  must  be  concluded,  therefore,  that  ideal  (100)  interfaces 
are  not  actually  formed.  It  is  possible  that  the  metallic  interface  band  is  eliminated  by 
interfacial  relaxation  or  reconstruction,  or  by  local  interdiffusion  and  reordering.21  It  is 
conceivable  that  some  macroscopic  polar  interfaces  are  composed  primarily  of  non-polar 
terraces  and  steps  on  an  atomic  scale. 

Recent  experimental  studies22*23  suggest  that  Ge/GaAs  interfaces  are  considerably  more 
complicated  chemically  and  structurally  than  was  originally  believed.  Because  of  the  similari¬ 
ty  in  atomic  size  of  Ge,  Ga,  and  As,  it  is  difficult  to  prevent  interdiffusion,  so  that  Interfaces 
tend  to  be  diffuse  rather  than  atomically  abrupt.  Moreover,  it  is  found  that  iaterfadal 
chemical  reactions  take  place,  and  that  the  interface  is  rough  on  an  atomic  scale.  Whan 
GaAs  is  grown  on  Ge,  antiphase  boundaries  can  be  formed  in  the  growing  GaAs,  further 
complicating  the  interfacial  structure.  Some  attempts  have  already  been  made  to  incorporate 
various  types  of  defects  into  interface  calculations.13  However,  much  remains  to  ha  done 
before  we  hare  an  adequate  understanding  of  the  effects  of  chemical  imparities,  vacancies, 
interstitials,  antisite  defects,  and  stoichiometric  mixing  on  interfacial  electronic  structure 
and  properties.  The  same  applies  for  misfit  dislocations,  stacking  disorder,  relaxation,  and 
reconstruction. 


III.  PASSIVATING  OXIDE  LAYERS:  Si/Si02  INTERFACES 

It  is  well  known  that  suitably  oxidized  silicon  substrates  are  protected  from  environmental 
contaminants  and  have  rather  low  densities  of  localized  interface  states  in  the  thermal  gap. 
Because  of  the  important  role  that  surface  passivation  plays  in  MOS  technology,24  there  to 
strong  motivation  for  understanding  the  physics  and  chemistry  of  the  Si/S102  interface.0"8 
In  this  section  we  will  be  concerned  with  the  atomic-scale  structure  of  the  Sl/Si02  inter¬ 
face,  including  the  residual  defects  that  occur  near  well-prepared  Si/Si02  interfaces  sad 
give  rise  to  residual  trapping  sites  and  fixed  charge. 

In  contrast  to  the  Ge/GaAs  interfaces  that  we  have  been  dealing  with  so  far,  the  Sl/Si02 
interface  connects  a  crystalline  substrate  (Si)  to  a  non-crystalline  overlayer  (vitreous  Si02, 
abbreviated  v-S102).  Whatever  the  nature  of  this  interface  actually  is  on  an  atomic  level,  the 
observed  low  density  of  interface  states  is  compelling  evidence  for  the  fact  that  nearly  all  of 
the  silicon  atoms  on  the  substrate  surface  have  their  bonds  saturated  one  way  or  another. 
How  does  this  come  about? 

Before  discussing  the  Si/Si02  interface  itself,  it  is  desirable  to  say  a  few  words  about  bulk 
v-S102.  It  is  generally  believed23  that  a  non-crystalline  material  such  as  v-Si02  can  be 
represented  by  a  continuous  random  network  (CRN)  of  S104  tetrahedra  joined  to  one 
another  at  their  common  oxygen  positions,  so  that  all  Si  and  O  bonds  are  saturated.  Com¬ 
puter  simulation  studies23,24  indicate  that  the  distributions  of  Si-O  bond  lengths  and  O-Si-O 
angles  are  relatively  narrow,  so  that  the  Si04  tetrahedra  maintain  their  structural  integrity 


Apptrilx  A  Semiconductor  Interfaces  18 

even  when  embedded  in  »-SIOj.  In  contrast,  the  distribution  of  Si-O-Si  bond  angles  involv- 
tng  adjacent  tetrahedra  is  relatively  broad.  Undoubtedly,  the  random  structure  of  the  Si04 
tetrabedra  in  v-SI02  is  made  possible  by  the  softness  of  the  Si-O-Sl  bond  angle.  The  fact 
that  the  SI04  tetrahedra  can  accommodate  themselves  to  different  SI02  crystal  structures 
mast  also  reflect  this  bond  angle  softness.  Ia  the  present  context,  this  softness  would  be 
expected  to  play  an  important  role  ia  accommodating  v-Si02  to  a  crystalline  Si  substrate. 

What  is  the  exact  nature  of  this  contact?  Is  the  silicon  substrate  atomically  smooth  or 
rough?  Does  it  have  steps  and  terraces?  Are  the  outermost  silicon  layers  strained  to 
accommodate  the  SK)2  overlayer?  Is  there  a  transition  region  between  the  Si  substrate  and 
the  "hulk”  T-SiOj  which  extends  over  several  atomic  layers,  or  is  the  transition  atomically 
abrupt?  If  there  is  a  gradual  change  la  average  composition,  described  by  the  stoichimetric 
f orarala  8K>X,  where  x  ranges  from  0  (at  Si)  to  2  (at  v-Si02),  what  are  the  bonding  ar¬ 
range  meats  that  contribute  to  this  average  composition?  We  do  not  have  the  space  to 
summarise  recent  experimental  evidence  hearing  on  these  questions,24  so  we  will  limit 
ourselves  to  dlarueslng  seme  ef  these  Issues  la  terms  of  two  recently  developed  theoretical 
models  which  are  highly  IdeaHsed  hat  nevertheless  instructive. 

The  first  of  these  is  a  batt-and -stick  model  developed  by  Pantelides  and  Long27  to  demon¬ 
strate  that  aa  atomically  abrupt  SI/SI02  Interfaces  can  be  formed  between  a  CRN  model  of 
v-SI02  and  the  (100)  face  of  Si.  Using  the  same  Si-O  bond  length  everywhere,  they  were 
able  to  connect  v-Si02  to  SI  on  the  scale  of  one  atomic  layer,  but  they  had  to  distort  many 
of  the  O-Si-O  and  Si-O-SI  bond  angles  near  the  Interface  considerably  to  avoid  introducing 
dangling  bonds.  The  successful  construction  of  this  ball-and-stick  model  dramatically 
dem  'rates  the  adaptability  of  v-SI02  to  a  silicon  substrate,  but  does  not  prove  that  such 
an  ate  really  abrupt  Si/SI02  Interface  would  actually  form  under  experimental  conditions. 

Since  the  more  distorted  bond  angles  are  energetically  unfavorable,  and  hence  unlikely  to  be 
realized  in  nature,  it  appears  necessary  to  Incorporate  two  additional  ingredients  into  the 
Pantelides-Long  model  in  order  to  make  it  more  realistic.  Ia  constructing  a  more  realistic 
version  of  this  model,  one  can  avoid  excessively  distorted  bond  angles  by  allowing  dangling 
bonds  to  appear,  and  one  can  saturated  these  dangling  bonds  by  O  and  H  atoms  and  OH 
groups,  as  undoubtedly  happens  during  actual  growth  conditions.28*29  The  extra  energy 
needed  to  support  dangling  bonds  would  be  largely  offset  by  the  saturation  of  these  bonds  by 
the  entities  just  mentioned.  (The  construction  of  a  model  embodying  these  features  is 
currently  being  investigated  at  our  laboratory  by  P.  Lambin  and  the  author.) 

The  net  effect  would  be  an  incomplete  connection  between  a  discontinuous  random  network 
of  v-Si02  and  the  silicon  substrate.  There  would  be  many  voids  of  varying  sizes  present  at  or 
near  the  interface,  but  most  of  the  dangling  bonds  defining  these  voids  would  be  saturated  by 
O  and  H  atoms  and  OH  groups.  Thus,  strain-relieving  voids  and  (OfH, OH) -saturated 
dangling  bonds  are  the  two  extra  ingredients  that  should  be  added  to  the  Panteiides-Long 
model  to  make  it  more  realistic.  The  local  stoichiometry  described  by  the  expression  SiOx 
would  of  course  be  determined  in  large  measure  by  the  O  atoms  and  OH  groups  which 
saturate  Si  bonds  and  represent  deviations  from  a  continuous  random  network  of  v-Si02. 

An  analogous  picture  emerges  from  a  totally  different  type  of  idealized  Si/Si02  interface 
model  that  we  developed  recently.28  In  this  other  model,  the  non-crystalline  Si02  overlayer 
is  replaced  by  an  idealized  crystalline  S102  overlayer  which  can  be  brought  into  exact 
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registry  with  the  Si  substrate.  If  one  begins  with  a  silicon  crystal*  expands  the  lattice  by  a 
factor  of  V2,  and  places  an  oxygen  atom  at  the  midpoint  of  each  dilated  Si-Si  bond,  one 
obtains  an  idealized  version  of  the  cristobalite  Si02  structure  having  very  nearly  the  same 
Si-O  bond  length  as  the  actual  structure,  but  straight  rather  than  crooked  Si-O-Si  bond 
angles.  This  is  a  reasonable  approximation  in  view  of  the  softness  of  the  Si-O-Si  bond 
angle.  One  can  now  obtain  exact  registry  between  the  (100)  face  of  Si  and  the  (100)  face  of 
idealized  cristobalite  SI02  by  rotating  one  of  these  45  degrees  with  respect  to  the  other  and 
bringing  them  together. 

The  Si  atoms  on  the  (100)  substrate  form  a  square  array  which  can  be  regarded  as  a  check¬ 
erboard  with  half  the  Si  atoms  occupying  "red"  squares  and  the  remaining  half  "black" 
squares.  At  the  idealized  Si/Si02  interface,  the  Si  atoms  occupying  the  "red"  squares  are 
shared  by  the  Si  substrate  and  the  Si02  overlayer,  while  the  Si  atoms  on  the  "black"  squares 
are  connected  only  to  the  Si  substrate,  having  two  dangling  bonds  each.  One  can  Improve  the 
realism  of  this  model  slightly  by  retaining  the  lattice  topology  and  allowing  the  atomic 
positions  to  relax,  but  this  does  not  change  the  essential  features  of  the  model. 

This  construction  leads  to  a  concrete,  easily  visualized  interface  model  which  emphasizes  that 
the  connection  between  the  Si  substrate  and  the  (idealized)  S102  overlayer  involves  only  half 
the  Si  substrate  atoms  directly.  The  remaining  half  of  the  Si  substrate  atoms  are  not 
attached  to  the  Si02  overlayer  at  alL  For  this  model  to  describe  an  actual  interface,  it  is 
again  essential  that  virtually  all  the  Si  substrate  dangling  bonds  are  saturated  by  H  or  O 
atoms  or  OH  groups  introduced  during  interface  formation.  Since  Si02  has  a  very  open 
structure,  it  is  reasonable  to  expect  such  entities  if  present  to  migrate  to  the  growing 
interface  and  saturate  the  dangling  bonds. 

Since  very  few  experiments  provide  direct  and  unambiguous  information  concerning  the 
exact  nature  of  the  atomic  structure  of  an  Si/Si02  interface,  it  is  essential  for  theoreticians 
to  develop  atomic-scale  models  that  can  be  used  as  a  guide  to  interpreting  experimental 
measurements.  The  two  models  already  discussed  are  clearly  idealizations,  but  they  are 
Instructive.  In  order  to  make  further  progress,  one  should  consider  the  atomic  processes  by 
which  the  interface  actually  forms,  rather  than  straightforward  attachments  between 
crystalline  Si0230  or  amorphous  S10227’31  to  a  Si  substrate. 

There  is  a  vast  literature  on  defects  in  Si02*~*  covering  interface  electronics  and  also 
radiation-damaged  glass.  It  is  encouraging  that  increasingly  refined  theoretical  models  have 
begun  to  appear  in  recent  years,32  as  well  as  experimental  studies  based  on  surface-sensitive 
spectroscopies.6-8,24  Nevertheless,  the  identification  of  the  defects  primarily  responsible  for 
trapping  and  recombination  at  Si/Si02  interfaces  is  still  largely  an  open  question.  Since  the 
thermal  (Si)  gap  is  only  1.1  eV  and  the  top  of  this  gap  lies  about  3.1  eV  below  the  top  of  the 
10  eV  forbidden  band  of  Si02,  most  of  the  defects  in  the  Si02  structure  would  be  expected 
to  give  rise  to  levels  that  He  outside  the  thermal  gap.  It  is  far  more  likely  that  defects  such 
as  dangling  Si  substrate  bonds  which  are  closely  coupled  to  the  Si  band  structure  give  rise  to 
levels  in  the  thermal  gap.  Support  for  this  view  comes  from  the  fact  that  the  deliberate 
introduction  of  hydrogen  or  deuterium  during  processing  removes  these  defect  levels, 
presumably  by  saturating  these  particular  dangling  bonds.28,2* 

Finally,  we  note  that  there  is  considerable  tailing  of  the  valence  and  conduction  bands  into 
the  thermal  gap.  This  tailing  is  supposedly  associated  with  the  disordered  atomic  structure  at 
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the  Si/Si(>2  interface,  but  the  details  are  still  not  understood.  Clearly,  a  great  deal  of 
additional  experimental  and  theoretical  work  will  have  to  be  done  before  we  have  a  more 
comprehensive  picture  of  this  interface. 


IV.  METAL-SEMICONDUCTOR  CONTACTS:  Pd2Si/Si(lll)  INTERFACES 

The  formation  of  ohmic  or  non-ohmic  metal-semiconductor  contacts  is  of  immense  techno¬ 
logical  importance.  During  the  past  few  years,  many  attempts  have  been  made  to  understand 
the  physical,  chemical,  metallurgical,  and  structural  factors  determining  Schottky  barrier 
heights  in  rectifying  contacts.9**0  Although  informative  chemical  trends  have  been  estab¬ 
lished  by  correlating  barrier  heights  with  electrochemical  data  for  many  systems,  the 
relationship  between  barrier  height  and  local  atomic  structure  at  the  interface  remains 
elusive.  Unfortunately,  the  atomic  arrangements  at  most  interfaces  are  too  poorly  charac¬ 
terized  to  provide  a  basis  for  serious  theoretical  analysis. 

For  most  metal-semiconductor  interfaces,  lattice-matched  constituents  and  atomically  abrupt 
interfaces  are  the  exception  rather  than  the  rule.  The  metal  and  semiconductor  regions  are 
usually  separated  by  a  disordered  transition  region  which  accommodates  the  discontinuity  in 
lattice  constant  as  well  as  the  transition  from  metallic  and  covalent  binding.  Most  of  these 
interfaces  are  spatially  diffuse  because  of  interdiffusion.  For  some  systems,  chemical 
reactions  may  take  place  between  the  metal  and  semiconductor,  leading  to  a  transition  region 
composed  of  ordered  or  disordered  compounds.  Such  non-ideal  systems  are  exceedingly 
difficult  to  treat  theoretically,  though  some  progress  can  be  made  by  considering  idealized 
lattice-matched  metal-semiconductor  interfaces.33 

Most  elementary  models  of  ideal  rectifying  contacts  are  based  on  early  ideas  by  Schottky  and 
Bardeen,  according  to  which  chemical  equilibrium  is  established  by  the  flow  of  electrons 
between  the  metal  and  the  semiconductor,  leading  to  bending  of  the  semiconductor  bands  and 
the  creation  of  a  rectifying  potential  barrier  at  the  interface.  The  absence  or  presence  of 
localized  interface  states  plays  a  crucial  role  in  these  models:  In  the  Schottky  limit 
(negligible  interface  states),  the  barrier  height  is  given  by  the  difference  between  the  metal 
work  function  and  the  semiconductor  electron  affinity.  In  the  Bardeen  limit,  there  are 
sufficiently  many  interface  states  present  to  pin  the  Fermi  level,  so  that  the  barrier  height  is 
determined  by  the  highest  occupied  interface  state,  and  is  largely  independent  of  the  metal 
work  function. 

If  the  Fermi  level  is  indeed  pinned  by  localized  interface  states  whose  energy  levels  lie  in  the 
forbidden  band  of  the  semiconductor,  how  do  these  interface  states  arise?  The  first  possibil¬ 
ity  is  that  they  are  similar  to  intrinsic  surface  states  on  the  semiconductor,  as  originally 
envisioned  by  Bardeen.  However,  in  many  semiconductors,  relaxation  and  reconstruction  at 
free  surfaces  push  surface  states  out  of  the  forbidden  band;33*34  the  same  could  happen  at 
metal-semiconductor  interfaces.33  Moreover,  Intrinsic  semiconductor  interface  states  can 
only  partially  resemble  intrinsic  surface  states  because  of  screening  by  the  adjacent  metallic 
region. 

A  second  possibility  is  that  electrons  can  be  localized  just  inside  the  semiconductor  by 
occupying  the  tails  of  the  wave  functions  which  spill  over  into  the  semiconductor  from  the 
metaL  These  are  known  as  metal-induced  gap  states,30  or  MIGS.  Such  localized  states  are 
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expected  to  occur  when  the  metal  is  at  least  several  atomic  layers  thick,  but  they  may  not 
occur  when  the  metal  is  only  one  atomic  layer  thick.  Since  the  Schottky  barrier  height  may 
already  have  its  limiting  ("bulk")  value  when  the  metal  coverage  is  monoatomic  or  sub- 
monoatomic  for  some  systems,  there  mav  be  no  direct  connection  between  MIGS  and 
Schottky  barrier  heights  for  such  systems.38 

A  third  possibility,  particularly  appropriate  for  "strong"  (reactive)  contacts,  is  that  the 
frontier  metal  and  semiconductor  atoms  form  chemical  bonds,37  and  that  the  energy  levels 
associated  with  these  chemical  bonds  or  hybridized  metal-semiconductor  states  extend  into 
the  forbidden  band.  A  fourth  possibility  is  that  the  localized  interface  states  arise  from 
structural  or  chemical  imperfections  in  the  semiconductor.38  Such  imperfections  include 
vacancies,  self-interstitials,  substitutional  and  interstitial  metal  impurities,  metallic  dusters, 
antisite  defects,  and  dislocations.  In  the  case  of  contacts  between  metals  and  hydrogenated 
amorphous  silicon,39  states  in  the  gap  already  present  because  of  structural  disorder  in  the 
silicon  are  expected  to  play  a  key  role  in  pinning  the  Fermi  level. 

Many  investigators  are  currently  making  a  concerted  effort  to  interpret  a  rapidly  increasing 
body  of  detailed  experimental  information  in  terms  of  different  defect  models.9*10*38,38  In 
view  of  the  fact  that  most  experimental  probes  do  not  have  the  spatial  resolution  necessary 
to  identify  structural  imperfections  unambiguously  on  an  atomic  scale,  a  general  consensus 
has  not  yet  been  reached  regarding  the  mechanisms  responsible  for  Schottky  barrier  forma¬ 
tion.  There  may  well  be  many  different  mechanisms  at  play  in  different  systems,  or  even  in 
the  same  system,  so  sorting  out  the  key  mechanisms  is  an  exceedingly  difficult  problem. 

The  demands  of  sub-micron  electronics  are  focusing  attention  on  rectifying  contacts  formed 
by  reacting  transition  metals  with  silicon.  There  is  considerable  interest  in  such  interfaces 
because  the  silicides  form  dose-packed  metallic  structures  which  are  impervious  to  bonding 
metals  such  as  AL  Thus,  diffuse  Al-Si  contacts  can  be  replaced  by  relatively  abrupt  Al- 
siiidde-Si  contacts.  Moreover,  some  of  the  silicides  form  epitaxial  layers  on  Si,  so  there  is  a 
greater  opportunity  to  synthesize  interfaces  which  are  more  nearly  ideal  and  hence  easier  to 
investigate  experimentally  and  theoretically.9*10,40 

There  have  already  been  a  number  of  experimental  and  theoretical  studies  of  the  electronic 
structure  of  bulk  silicides.40  Some  attempts  are  currently  being  made  to  study  the  electronic 
structure  of  silicide-silicon  interfaces  theoretically,  taking  advantage  of  their  epitaxial 
character.41*42  In  the  following,  we  will  discuss  the  essential  features  of  a  model  we  have 
recently  developed  for  understanding  Schottky  terriers  at  Pd2Si/Si(lll)  interfaces.42  This 
model  is  of  general  interest  because  it  goes  beyond  the  idea  of  localized  interface  states 
being  solely  responsible  for  Schottky  terrier  heights. 

According  to  our  model,  Pd  impurities  diffuse  into  the  Si  substrate  in  advance  of  the  growing 
interface,  forming  triangular  clusters  between  pairs  of  Si  vacancies.43  Since  these  clusters 
are  more  easily  accommodated  in  hexagonal  (2H)  Si  than  in  cubic  (3C)  Si,  the  presence  of 
these  clusters  induces  a  hexagonal  Si  transition  region  between  the  cubic  SI  substrate  and  the 
Pd2Si  overlayer.  The  measured  Schottky  terrier  height  (0.71  to  0.73  eV)  is  readily  account¬ 
ed  for  by  noting  first  that  the  forbidden  tend  width  of  the  hexagonal  Si  transition  region  is 
about  0.85  eV,  and  second  that  the  Fermi  level  is  pinned  slightly  above  the  Si  valence  tend 
edge  by  hybridized  Pd-Si  states  at  the  interface.41*42  The  same  overall  effect  could  be 
achieved  by  having  stacking  disorder  in  the  Si  near  the  interface,  rather  than  a  coherent 
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hexagonal  crystal  structure.  We  would  again  argue  that  the  stacking  disorder  is  induced  by 
the  presence  of  the  triangular  Pd  impurities. 

In  short,  our  Pd2Si/Si(lll)  model  suggests  that  Schottky  barrier  heights  can  be  determined 
not  only  by  interface  states  which  pin  the  Fermi  level,  but  also  by  the  presence  of  a  different 
crystal  structure  or  stacking  disorder  near  the  interface.  Such  structural  features  — 
coherent  or  random  —  can  affect  the  host  band  structure-  modifying  the  effective  semicon¬ 
ductor  band  gap  at  the  interface,  and  hence  the  Schottky  barrier  height  as  well.  Incidentally, 
stacking  faults  have  also  been  invoked  recently  to  account  for  certain  reconstructions  of  free 
silicon  surfaces.44 


V.  CONSTRUCTION  OF  ATOMIC-SCALE  INTERFACE  MODELS 

One  popular  theoretical  approach  for  determining  the  atomic  structure  of  interfaces  (and 
surfaces)  involves  calculating  the  electronic  structure  for  a  number  of  plausible  models,  and 
then  determining  which  set  of  results  agrees  most  closely  with  experiment.  *3,34  Even 
assuming  that  the  experimental  data  are  reproducible  and  correctly  interpreted,  this  approach 
may  overlook  physically  significant  atomic  arrangements  altogether,  and  in  the  end  not 
provide  insight  as  to  why  one  particular  model  gives  a  better  account  of  experiment  than 
another. 

Another  popular  theoretical  approach  involves  calculating  the  total  energy  of  the  system,  and 
then  determining  the  equilibrium  geometry  by  minimizing  the  total  energy  as  a  function  of 
atomic  positions.  This  is  a  very  powerful  approach,  and  has  been  applied  with  great  success 
recently  to  a  number  of  relatively  simple  systems.13,34  Because  of  the  extensive  computa¬ 
tional  effort  involved,  it  is  unlikely  that  this  approach  will  soon  encompass  the  more  compli¬ 
cated  systems  that  are  of  primary  experimental  interest. 

So  it  is  essential  to  look  for  still  other  approaches,  bearing  in  mind  that  further  progress  in 
determining  the  electronic  structure  of  interfaces  is  severely  limited  by  our  ability  to 
construct  realistic  atomic  models  of  interfaces.  Instead  of  assuming  plausible  atomic 
positions  at  the  outset  and  then  relaxing  these  positions  so  as  to  minimize  the  total  energy, 
keeping  the  topology  or  connectivity  the  same,  we  believe  theoreticians  should  focus  directly 
on  the  atomic  processes  by  which  interfaces  are  formed  by  the  progressive  accumulation  of 
atoms  and  molecules.  We  should  attempt  to  simulate  the  growth  of  covalently  bonded 
crystals  at  the  atomic  level  using  computers,  perhaps  interactively.  The  challenge  is  to 
establish  algorithms  that  will  guide  this  growth  along  physically  and  chemically  reasonable 
paths,  possibly  taking  advantage  of  concepts  arising  in  the  field  of  artificial  intelligence. 

Computer  simulation  studies43  have  already  provided  us  with  considerable  insight  into  the 
nature  of  crystal  growth,  but  such  studies  have  been  confined  almost  exclusively  to  systems 
such  as  metals  and  simple  liquids  which  can  be  described  by  hard  sphere  interactions.43 
Current  attempts  at  modeling  covalently  bonded  systems  are  still  rather  primitive,47,48  but 
this  avenue  of  theoretical  research  is  of  the  greatest  importance,  and  deserves  considerably 
more  attention  than  it  has  received  so  far. 
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VL  FORMATION  OF  NOVEL  INTERFACES  AND  SUPERLATTICES 

During  the  next  few  years  we  anticipate  increasing  experimental,  theoretical,  and  computa¬ 
tional  efforts  directed  at  understanding  the  growth  and  stabilization  of  novel  interfaces  and 
superlattices,  with  particular  emphasis  on  the  underlying  atomic  processes: 

•  Amorphous  overlayers  on  crystalline  substrates.  —  We  have  already  discussed  the  Si/Si02 
interface  as  an  important  example  of  a  crystalline-amorphous  interface.  A  somewhat  simpler 
example  is  the  interface  between  crystalline  Si  and  amorphous  Si.  Recent  activity  in  laser 
processing  and  laser  annealing49  has  focused  attention  on  the  fundamental  processes 
associated  with  localized  melting  and  subsequent  solidification  of  Si  into  the  amorphous 
phase,  as  well  as  with  recrystallization  of  amorphous  Si.  Some  attempts  have  already  been 
made  to  model  crystalline-amorphous  Si  transformations  as  well  as  the  growth  of  Si02  on  Si 
on  an  atomic  scale,50  but  progress  thus  far  has  been  minimal.  We  don’t  really  know  very 
much  about  the  dynamical  processes  responsible  for  the  interconversion  between  crystalline 
and  amorphous  phases  of  covalently  bonded  solids.  If  we  could  understand  the  atomic 
dynamics  by  which  crystals  transform  into  amorphous  structures  and  vice  versa,  we  would  be 
able  to  determine  the  atomic  arrangements  at  crystalline-amorphous  interfaces.  This  is 
clearly  a  challenge  for  the  future. 

•  Metastable  overlayers  on  solid  substrates.51  —  Perhaps  the  best-known  examples  are 
strained-layer  interfaces  and  superlattices.11  These  metastable  structures  have  many 
attractive  physical  properties,  so  they  are  being  intensively  studied.  An  important  challenge 
is  to  find  ways  of  stabilizing  such  structures  so  that  they  can  perform  electronic  functions 
under  widely  varying  conditions.  Other  examples  are  low-temperature  phase  a  tin  ("grey 
tin")  and  a-Sn:Ge  alloys  grown  on  InSb  and  CdTe  substrates.52  Still  other  examples  are 
single-crystal  (GaAs)|,xGe2x  alloys,  whose  band  gaps  vary  non-monotonicaliy  with  compos¬ 
ition  in  a  most  remarkable  manner.55  This  behavior  has  been  interpreted  as  evidence  for  an 
order-disorder  transition.54 

•  Polytypism,55  stacking  disorder,56  and  random  superstructures.57  —  The  boundary 
between  cubic  ZnS  (sphalerite)  and  hexagonal  ZnS  (wurtzite)  is  a  classic  example  of  a 
polymorphic  interface.58,59  Long-period  polytypes  of  SiC  and  ZnS  can  be  rightfully  regarded 
as  structurally  modulated  superlattices.  The  mechanisms  leading  to  the  growth  of  such 
superlattices  are  still  not  understood,  in  spite  of  decades  of  study.  Admixed  random  and 
non-random  stacking  sequences  of  ZnS  can  be  formed  by  evaporating  thin  films  of  ZnS  on 
suitable  substrates.  These  films  are  composed  of  random  striations  of  cubic  and  hexagonal 
ZnS.58,60  Such  films  are  of  interest  because  they  can  be  used  to  generate  photovoltages 
considerably  larger  than  the  band  gap  of  ZnS.  Another  challenge  for  the  future  is  learning 
to  synthesize  structurally  modulated  heterostructures  as  well  as  we  have  already  learned  to 
synthesize  compositionally  modulated  heterostructures1’2,5,11  and  modulation  doped 
heterostructures.5  Hopefully,  an  improved  theoretical  understanding  of  interface  formation 
will  accelerate  progress. 
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VII.  CONCLUDING  REMARKS 

Recent  experimental  studies  have  provided  a  great  deal  of  important  information  concerning 
the  nature  of  semiconductor  interfaces.  Many  experimental  results  can  be  readily  interpret* 
ed  in  terms  of  electronic  structure  pictures  based  on  idealized  interface  geometries.  Howev¬ 
er,  it  has  proved  difficult  to  elucidate  the  atomic  arrangements  at  interfaces  unambiguously 
because  even  the  most  refined  surface-sensitive  spectroscopies  lack  the  necessary  resolution. 
Further  theoretical  progress  on  electronic  structure  and  more  incisive  interpretation  of 
experimental  measurements  both  await  a  better  understanding  of  the  atomic  structure  of  real 
interfaces.  The  development  of  realistic  atomic  structure  models  is  a  challenging  theoretical 
problem  which  is  ideally  suited  for  imaginative  computational  analysis. 
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ELECTRONIC  STRUCTURE  OF  THE  Si/S102  INTERFACE 

Frank  Herman 
IBM  Research  Laboratory 
San  Jose,  CA  95193  USA1 

ABSTRACT.  After  sketching  various  theoretical  models  that  have  been  proposed  for 
describing  the  Si/Si02  interface  on  an  atomic  scale,  we  will  discuss  an  idealized  model  which 
is  designed  to  simulate  the  average  contact  between  a  silicon  crystal  and  its  oxide  overlayer. 
According  to  this  model,  the  Si/S102  interface  is  represented  by  the  boundary  between  two 
crystalline  domains,  the  first  being  the  silicon  substrate,  and  the  second  an  idealized  crystal¬ 
line  form  of  Si02,  diamond-like  beta  cristobalite.  Nearly  perfect  registry  between  Si  and 
Si02  is  obtained  by  placing  the  (100)  face  of  the  former  next  to  the  45°  rotated  (100)  face 
of  the  latter.  Half  the  Si  atoms  at  the  interface  are  four-fold  coordinated,  connecting  the  Si 
and  Si02  regions  to  one  another.  The  remaining  half  are  connected  only  to  the  silicon 
substrate  and  are  thus  two-fold  coordinated.  Localized  electronic  states  in  the  thermal  gap 
at  the  interface  are  associated  with  the  dangling  bonds  at  the  two-fold  coordinated  Si  atoms. 
If  H,  O,  or  OH  groups  are  attached  to  these  Si  atoms,  eliminating  the  dangling  bonds,  the 
thermal  gap  interface  states  are  also  eliminated.  Even  though  this  model  ignores  structural 
disorder,  it  does  provide  considerable  insight  into  the  nature  of  the  Si/Si02  interface.  We 
will  discuss  some  of  the  implications  of  this  model,  as  well  as  the  results  of  model-based 
theoretical  calculations.  In  particular,  we  will  consider  the  electronic  structure  of  the 
interface  including  the  effects  produced  by  the  introduction  of  Si  and  O  vacancies  and 
interstitials  and  other  types  of  localized  defects. 


1.  INTRODUCTION 

Because  of  the  important  role  played  by  the  Si/S102  interface  in  modern  semiconductor 
technology,  considerable  experimental  and  theoretical  effort  has  been  devoted  to  the  study  of 
this  interface.2-8  Oxidation  of  a  silicon  surface  produces  an  amorphous  oxide  layer  whose 
atomic  arrangement  is  not  well  understood.  Many  authors  have  suggested  that  there  is  a 
non-stoichiometric  transition  layer,  often  described  as  SiOx,  separating  the  silicon  substrate 
from  the  stoichiometric  Si02  layer.  (Here  x  lies  between  1  and  2.)  The  detailed  atomic 
structure  of  the  transition  layer,  including  its  uniformity  and  width,  its  structural  disorder,  its 
affinity  for  impurity  atoms,  etc.,  continues  to  be  a  subject  of  debate.  Of  course,  the  nature 
of  the  transition  layer  is  expected  to  depend  on  the  manner  in  which  the  oxide  layer  is 
formed.  By  the  same  token,  a  better  understanding  of  the  transition  layer  could  lead  to 
further  improvements  in  the  preparation  of  Si/S102  interfaces. 

The  technological  importance  of  the  Si/S102  interface  arises  from  the  fortunate  circum¬ 
stance  that  suitably  oxidized  silicon  substrates  are  protected  from  environmental  contami¬ 
nants  and  have  rather  low  densities  of  localized  Interface  states  in  the  thermal  gap.  There 
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continues  to  be  considerable  interest  in  the  nature  and  identity  of  the  residual  defects  that 
occur  near  well-prepared  interfaces,  giving  rise  to  residual  trapping  sites  and  fixed  charge. 
Many  imaginative  models  of  defects  at  the  SI/SI02  interface  and  in  Si02  have  been 
proposed,2-11  but  relatively  few  of  these  defects  have  been  unambiguously  identified,  even 
by  such  precise  measurements  as  ESR.11  In  view  of  the  difficulties  encountered  in  practice 
in  attempting  to  interpret  most  experimental  measurements  in  terms  of  specific  defects,  we 
thought  it  would  be  useful  to  develop  a  simple  model  for  the  Si/Sl02  interface  which  could 
provide  a  conceptual  framework  for  thinking  about  interface  defects  as  well  as  a  theoretical 
framework  for  carrying  out  numerical  calculations. 

In  the  first  paper  on  this  subject,12  we  described  the  construction  of  idealized  Sl/S102 
interfaces  and  superlattices  and  demonstrated  by  first  principles  band  structure  calculations 
that  there  are  no  localized  interface  states  in  the  thermal  gap  if  all  the  Si  atoms  at  the 
interface  are  saturated.  In  a  second  paper,12  we  showed  how  localized  defects  in  S102  could 
be  Investigated  theoretically  by  introducing  periodic  arrays  of  non-interacting  defects  in  a 
Si02  supercell  normally  containing  24  atoms.  In  a  third  paper,14  we  introduced  periodic 
arrays  of  non-interacting  defects  in  idealized  Si/Si02  superlattices,  thereby  simulating  the 
presence  of  localized  defects  at  or  near  the  Si/Si02  interface.  In  the  present  paper  we  will 
first  mention  some  recent  attempts  by  others  to  develop  a  theoretical  model  of  S102  and  the 
Si/Si02  interface,  and  then  we  will  discuss  our  own  model  briefly  and  summarize  some  of 
our  principal  conclusions. 


2.  THEORETICAL  MODELS  OF  THE  Si/Si02  INTERFACE 

The  extensive  literature  on  crystalline  and  non-crystalline  Si02  and  on  Si/Si02 
interfaces2-11  contains  a  variety  of  models  of  the  oxide  layer  and  the  Si/Si02  interface.  In 
this  paper,  we  will  focus  on  atomic  scale  rather  than  on  phenomenological  models.  Nearly  all 
the  models  that  have  been  proposed  for  describing  the  atomic  arrangements  at  the  interface 
are  speculative.  Some  of  the  proposals  are  quite  specific.  They  assume  particular  types  of 
linkages  between  Si  and  Si02,  well-defined  structural  imperfections  arising  from  modifica¬ 
tions  in  the  normal  linkages,  as  well  as  chemical  imperfections  due  to  the  presence  of  specific 
foreign  species.  The  remaining  proposals  tend  to  be  descriptive  or  impressionistic,  invoking 
dislocations,  disjunctions,  varying  degrees  of  stoichiometry  and  crystallinity,  etc.  The  first 
set  of  proposals  can  be  distinguished  from  the  second  by  the  ease  with  which  ball-and-stick 
models  can  be  constructed  to  illustrate  the  proposed  arrangements. 

There  has  also  been  considerable  debate  regarding  the  initial  formation  of  the  oxide  layer, 
and  how  this  layer  might  change  as  it  grows.2-8  Some  authors  believe  that  a  better  under¬ 
standing  of  the  initial  stages  of  oxidation  could  provide  valuable  clues  concerning  the  nature 
of  the  ultimate  oxide  layer;  other  authors  disagree.  In  any  event,  it  is  very  important  to  have 
heuristic  atomic  models  of  the  formation  and  ultimate  structure  of  the  Si/Si02  interface, 
because  such  ideas  could  lead  to  improved  processing  techniques,  new  device  concepts,  and 
better  understanding  of  still  more  complex  semiconductor/insuiator  interfaces. 

Up  to  now,  most  Si/Si02  interface  modeling  studies  have  been  qualitative,  being  based  on 
simple  physical  and  chemical  reasoning.  Many  of  the  arguments  put  forward  In  favor  of 
specific  atomic  models  seem  qnite  plausible,  but  it  is  usually  difficult  to  verify  these  interface 
models  experimentally.  Therefore,  a  constructive  step  would  be  to  study  some  of  these 
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models  theoretically  and  obtain  quantitative  estimates  of  crucial  electronic  properties.  If  one 
wishes  to  carry  out  detailed  electronic  structure  calculations  for  model  systems,  it  is  obvious* 
ly  necessary  to  have  specific  atomic  arrangements  in  mind.  The  essential  difficulty  is  that  we 
do  not  know  the  actual  structure  of  the  Si/S102  interface  on  an  atomic  scale. 

The  overall  problem  can  be  dealt  with  in  three  parts:  In  the  first  and  easiest  part,  one  makes 
certain  assumptions  about  the  silicon  substrate.  For  example,  one  chooses  a  specific 
crystalline  orientation,  and  then  assumes  that  the  silicon  is  structurally  perfect  up  to  the 
interface.  Relaxation  and  reconstruction  effects  can  be  deferred  to  a  later  stage  or  express¬ 
ly  ignored.  The  second  part  concerns  the  description  of  the  oxide  film  proper,  allowing  if 
possible  for  its  non-crystalline  nature.  In  dealing  with  non-crystalline  Si02,  one  can  build  on 
our  improved  knowledge  of  amorphous  and  glassy  semiconductors.9’15"19  The  third  part,  the 
most  difficult  of  all,  involves  connecting  the  silicon  substrate  to  the  non-crystalline  oxide 
layer  in  a  physically  and  chemically  plausible  manner. 

One  popular  model  for  non-crystalline  Si02  and  the  Si/Si02  interface  is  the  cluster  Bethe 
lattice  model,20’21  which  takes  topological  aspects  of  structural  disorder  into  account  in  an 
elegant  manner.  Because  this  model  is  rather  simple  in  form  and  can  be  easily  parameter¬ 
ized,  it  can  be  used  to  study  electronic  and  lattice  vibrational  spectral  properties  associated 
with  ideal  and  defective  interfaces.  But  because  these  models  are  based  on  adjustable 
parameters  that  are  assigned  ad  hoc  values,  they  are  not  particularly  trustworthy  from  a 
quantitative  point  of  view.22  Nevertheless,  they  do  have  considerable  pedagogical  value. 

Another  nopular  model  of  non-crystalline  S102  is  the  continuous  random  network 
model,15"19  according  to  which  Si02  is  composed  of  a  randomly  oriented  collection  of  S104 
tetrahedra  which  are  linked  to  one  another  to  form  a  continuous  network  of  Si-O-Si  bonds. 
It  is  assumed  that  the  SIO4  tetrahedra  tend  to  maintain  their  structural  integrity  (bond 
lengths  and  bond  angles)  in  different  environments,  and  that  the  necessary  structural 
versatility  is  provided  by  the  Si-O-Si  bond,  whose  angle  may  assume  a  wide  range  of  values 
at  little  cost  in  energy.4  Models  of  this  type  have  recently  been  used  to  study  bulk  Si0222  as 
well  as  the  Si/Si02  interface.24  In  the  study  of  bulk  Si02,  cyclic  boundary  conditions  were 
imposed  on  a  large  disordered  molecular  S102  cluster,  producing  a  computationally  tractable 
model  of  non-crystalline  Si02.  The  object  here  was  to  investigate  electron  and  hole 
localization  induced  by  the  structural  disorder. 

In  the  study  of  the  Si/Si02  interface,  it  was  demonstrated  that  a  rather  abrupt  interface 
could  be  constructed  within  the  framework  of  a  continuous  random  network  model,  a  rather 
important  result.  Unfortunately,  the  resulting  atomic  model  is  much  too  complex  for 
detailed  electronic  structure  studies.  In  both  Refs.  23  and  24,  the  construction  of  continuous 
random  network  models  is  based  on  ideas  regarding  the  nature  of  non-crystallinity  which  are 
highly  subjective  and  which  may  not  necessarily  describe  the  experimental  situation.  It  is 
possible  that  non-crystalline  Si02  near  a  Si/Si02  interface  does  not  actually  form  a  com¬ 
pletely  continuous  random  network.  That  is  to  say,  the  network  may  contain  discontinuous 
elements  such  as  non-bridging  O  atoms  or  three-fold  coordinated  Si  atoms  connected  to  H 
or  OH  groups,  or  voids  composed  of  aggregates  of  such  discontinuous  elements.19 

There  are  many  theoretical  efforts  currently  underway18’19’25  whose  aim  is  to  formulate 
general  guiding  principles  for  constructing  atomic  models  of  non-crystalline  materials. 
Combined  with  intuitive  notions20  and  computer  modeling,27*28  these  theoretical  studies  may 
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eventually  lead  us  to  objective  atomic  models  which  we  can  use  with  confidence  for  electron¬ 
ic  structure  studies  of  non-crystalline  Si02.  In  the  meanwhile,  the  atomic  architecture  of  the 
Sl/S102  interface  remains  an  open  question. 

We  already  know  from  extensive  experimental  and  theoretical  studies2,3*29  that  the  gross 
features  of  the  electronic  structure  of  non-  crystalline  and  various  crystalline  forms  of  Si02 
are  essentially  the  same,  reflecting  the  common  Si04  building  block.  Only  the  finer  details 
of  electronic  structure  are  affected  by  the  special  manner  in  which  the  Si04  tetrahedra  are 
linked  together  in  a  particular  structural  modification.  In  view  of  the  overwhelming  impor¬ 
tance  of  the  short  range  order  (linked  Si04  tetrahedra)  in  determining  the  essential  features 
of  the  electronic  structure  of  Si02,  and  the  likelihood  that  the  long  range  disorder  affects 
only  the  finer  details,  we  will  explicitly  ignore  the  structural  disorder  of  Si02  and  concen¬ 
trate  on  models  of  the  interface  which  represent  the  average  contact  between  Si  and  Si02. 


3.  IDEALIZED  Si/Si02  SUPERLATTICES  AND  INTERFACES 

Our  overall  program  is  as  follows.  First,  we  will  disregard  structural  disorder  in  Si02,  and 
represent  S102  by  a  suitably  chosen  periodic  structure.  We  will  then  form  ordered  Si/Si02 
interfaces  by  attaching  crystalline  Si02  to  crystalline  Si.  Next,  we  will  construct  a  Si/Si02 
superlattice  by  stacking  Si  and  Si02  slabs  on  top  of  one  another.  The  slabs  will  be  made 
sufficiently  thick  so  that  successive  interfaces  will  not  interact  with  one  another.  Because 
the  Si/S102  superlattice  is  a  periodic  structure  in  three  dimensions,  we  can  use  the  highly 
developed  methods  of  band  theory30  to  calculate  the  electronic  structure.  The  results  of  such 
calculations  include  a  description  of  localized  interface  states,  including  their  energy  levels 
and  charge  distributions.  Finally,  we  can  introduce  structural  and  chemical  imperfections, 
and  study  their  effects  on  the  electronic  structure. 

To  carry  out  the  above  program,  we  begin  by  representing  the  oxide  layer  by  an  idealized 
crystalline  form  of  Si02  having  as  simple  a  unit  cell  as  possible.  For  this  purpose  we  will  use 
beta  cristobaiite,31  except  that  we  will  straighten  out  the  Si-O-Si  bonds  in  the  actual 
structure,  making  them  ail  linear  rather  than  crooked  (all  bond  angles  180°  rather  than  about 
140°).  Assuming  an  Si-O  bond  length  of  1.633  A,32  the  linear  Sl-O-Si  bond  length 
becomes  twice  this,  and  the  unit  cube  edge  of  idealized  diamond-like  Si02  becomes  7.543  A, 
which  is  5  percent  larger  than  the  unit  cube  edge  of  actual  beta  cristobaiite,  7.16  A.  This  5 
percent  difference  reflects  the  contraction  produced  by  the  bent  bonds  in  the  actual  struc¬ 
ture.  Except  as  otherwise  noted,  we  will  henceforth  regard  Si02  as  having  the  idealized 
diamond-like  beta  cristobaiite  structure  with  linear  Si-O-Si  bonds. 

It  is  an  interesting  exercise  to  find  ways  of  attaching  crystalline  Si02  to  crystalline  Si  so  that 
as  many  bonds  as  possible  are  formed.  The  simplest  arrangement  of  all  is  to  place  the  (100) 
face  of  Si  next  to  the  45°-rotated  face  of  (100)  Si02.  This  produces  almost  perfect  registry 
because  the  ratio  of  the  lattice  constants  of  Si02  and  Si  is  very  nearly  equal  to  the  square 
root  of  2  (7.543/5.431  »  1.39).  To  put  it  more  directly,  we  obtain  nearly  perfect  registry 
in  this  fashion  because  the  ratio  of  the  straight  Si-O-Si  bond  length  in  Si02  to  the  Si-Si 
bond  length  in  silicon  is  very  nearly  equal  to  the  square  root  of  2.  (For  actual  calculations, 
we  will  increase  the  Si-O-Si  bond  length  slightly  so  that  the  registry  is  exact.) 
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We  thus  obtain  an  extremely  simple  model  of  the  Si/Si02  interface,  namely,  a  fully  coher¬ 
ent,  ideally  abrupt  Interface  between  silicon  and  crystalline  S102.12"14  This  construction 
leads  to  an  interface  at  which  half  the  Si  substrate  atoms  are  four-fold  coordinated,  the 
remaining  half  being  two-fold  coordinated,  with  two  dangling  bonds  each.  We  can  describe 
the  interface  as  a  checkerboard  with  Si  atoms  common  to  Si  and  Si02  occupying  the  red 
squares,  and  unsaturated  Si  atoms  (two  dangling  bonds  each)  occupying  the  black  squares.  It 
is  possible  to  use  other  Si  and  Si02  crystal  faces  and  join  them  together  in  a  similar  manner, 
but  these  alternate  constructions  lead  to  more  complex  Interface  matching  patterns, 
stretched  Si-O-Si  bonds,  and  larger  repeat  periods,  making  them  less  suitable  for  modeling. 
It  is  Interesting,  nevertheless,  that  even  these  more  complex  models  exhibit  the  property  that 
roughly  half  the  Si  substrate  atoms  are  attached  to  the  Si02,  while  the  remaining  half  are 
unattached.  If  we  are  to  adopt  an  idealized  model,  we  might  as  well  adopt  the  simplest  one 
possible,  and  this  is  unquestionably  the  (100)  4S°-rotated  geometry  already  described. 


Because  of  the  43°  rotation  at  the  Interface,  the  four  bonds  emanating  from  the  common  Si 
atoms  deviate  from  ideal  tetrahedral  geometry.  Earlier  studies13  suggest  that  the  silicon 
lattice  is  considerably  more  rigid  than  the  S102  lattice,  so  that  most  of  the  strain  would  be 
relieved  by  lattice  relaxation  in  the  Si02  region.  We  have  made  some  progress  in  taking 
lattice  relaxation  effects  into  account,  but  on  the  basis  of  our  work  to  date13  we  view  the 
distortion  from  ideal  tetrahedral  geometry  at  the  common  Si  atoms  as  a  relatively  minor 
shortcoming  of  our  model.  The  thermal  (Si)  band  gap  is  determined  in  part  by  the  common 
Si  atoms  and  the  fact  that  they  are  saturated.  The  deviation  from  local  tetrahedral  geometry 
appears  to  be  only  of  secondary  importance. 


It  is  possible  to  avoid  the  lattice  strain  associated  with  the  45°  rotation  altogether  by 
considering  a  more  exotic  interface  model.  Here  each  common  Si  atom  is  replaced  by  a 
Si -Si  molecule,  one  end  of  which  terminates  the  silicon  substrate,  and  the  other  end  the 
SI02  region.  The  45°  rotation  is  then  absorbed  by  the  Si -Si  double  bonds.  The  net  result 
is  that  all  bond  angles  at  the  interface  and  elsewhere  are  tetrahedral.  In  this  model  the  Si-Si 
bond  length  is  2.352  A  in  the  silicon  crystal  and  slightly  less  (2.252  A)  in  Si-Si32.  We  have 
not  adopted  this  more  exotic  model  for  two  reasons.  First,  the  Si-Si  bond  is  a  weak  one,33 
so  there  is  some  question  as  to  whether  it  would  be  stable  in  the  environment  just  described. 
Secondly,  we  have  carried  out  comparison  calculations  using  the  more  exotic  as  well  as  the 
original  model,  and  we  find  that  the  Si  double  bonds  lead  to  a  significant  reduction  in  the 
Si/Si02  thermal  gap,  which  is  contrary  to  experiment.  So  we  will  stick  with  the  original 
model. 


We  do  not  regard  the  use  of  straight  rather  than  crooked  Si-O-Si  bonds  in  Si02  as  a  serious 
defect  in  our  modeL  We  have  already  said  that  there  is  very  little  energy  required  to  change 
the  Si-O-Si  bond  angle  from  140  to  180°. 4  Some  bond  distortion  is  inevitable  at  the 
Si/S102  interface  anyway.  Moreover,  the  Si-O-Si  bond  angles  are  not  expected  to  play  a 
major  role  in  determining  the  electronic  structure.29 

By  using  a  periodic  model  and  introducing  an  idealized  interface,  we  have  focused  attention 
on  the  average  contact  between  Si  and  S102,  highlighting  the  fact  that  this  average  contact 
requires  approximately  half  the  Si  substrate  atoms  to  be  attached  to  the  Si02,  and  the 
remaining  half  to  be  unattached.  The  unattached  (two-fold  coordinated)  Si  atoms  face  onto 
cavities  in  the  adjoining  Si02  region,  presumably  making  it  easy  for  species  such  as  O,  H,  or 
OH  to  find  their  way  to  the  beckoning  dangling  bonds.  Because  of  the  periodicity,  and  the 
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simple  manner  in  which  we  have  been  able  to  attach  SI02  to  the  silicon  substrate,  we  can 
describe  various  forms  of  the  transition  layer,  so-called  non-stoichiometric  SiOx,  on  an 
atomic  scale.  According  to  our  model,  this  transition  layer  does  not  actually  exist  if  each 
two-fold  coordinated  Si  atom  at  the  interface  is  saturated  by  two  H  atoms.  That  is  to  say, 
the  passage  from  Si  to  Si02  is  essentially  abrupt.  By  saturating  some  or  all  of  these  two-fold 
coordinated  Si  atoms  with  O  or  OH  groups,  we  can  introduce  varying  O  concentrations  in 
one  atomic  layer  just  beyond  the  silicon  substrate,  this  now  becoming  the  transition  layer. 
Perhaps  simple  geometrical  pictures  such  as  these  will  prove  more  helpful  in  interpreting 
experimental  measurements  than  complex  models  of  a  non-stoichiometric  transition  layer. 


4.  ELECTRONIC  STRUCTURE  CALCULATIONS  FOR  INTRINSIC  INTERFACES 

Having  constructed  an  idealized  interface  model,  we  proceed  to  study  its  electronic  structure 
by  considering  an  Si/Si02  superlattice  composed  of  alternating  Si  and  Si02  slabs  which  are 
sufficiently  thick  to  isolate  adjacent  interfaces  from  one  another.  Most  of  our  studies  are 
based  on  a  21  atom  per  unit  cell  superlattice  consisting  of  5  Si  layers  (2  atoms  per  layer) 
alternating  with  3  tiers  of  Si04  tetrahedra  (II  atoms).  This  superlattice  with  two  additional 
O  atoms  forming  Si**0  double  bonds  with  the  two-fold  coordinated  Si  atoms  is  shown  in.*2 

Our  studies  are  based  on  the  first-principles  extended  muffin  tin  orbital  (EMTO)  method.34 
All  of  the  calculations  reported  in  Refs.  12  to  14  and  summarized  below  are  based  on 
superlattice  charge  distributions  constructed  from  neutral  free  atom  self-consistent  charge 
densities.33  Although  the  superlattice  charge  distributions  themselves  are  not  exactly 
self-consistent,  such  calculations  give  a  good  account  of  Si,  Si02,  and  many  other  types  of 
semiconductors  and  insulators.  A  fully  self-consistent  EMTO  computer  program  has 
recently  been  developed,36  and  this  program  will  be  used  in  subsequent  studies  of  Si/Si02 
interfaces.  We  believe  that  the  EMTO  calculations  done  to  date  on  Si/Si02  provide 
reasonable  estimates  of  the  electronic  structure  of  ideal  and  defective  Sl/Si02.  Refinements 
based  on  fully  self-consistent  EMTO  calculations  will  be  reported  in  due  course. 

The  essential  results  for  an  intrinsic  Si/S102  superlattice  can  be  summarized  as  follows:  (a) 
Localized  interface  states  occur  in  the  thermal  gap  of  Si  if  the  dangling  bonds  are  left 
unsaturated,  (b)  These  localized  states  are  removed  from  the  thermal  gap  if  the  dangling 
bonds  are  saturated  by  pairs  of  H  atoms,  single  O  atoms,  OH  groups,  etc.  These  conclusions 
apply  to  unrelaxed  lattice  geometries.  We  have  made  some  attempts  to  relieve  the  lattice 
strain  using  Monte  Carlo  methods.*3.  Because  of  the  limited  width  of  the  Si02  region  and 
the  severe  constraints  imposed  by  the  rigid  intervening  Si  regions,  the  Si02  regions  relax 
only  slightly,  leading  to  negligible  changes  in  the  energy  level  structure  and  to  no  change  in 
conclusions  (a)  and  (b)  above. 


5.  Si  AND  O  VACANCIES  AND  INTERSTITIALS 

Within  the  framework  of  a  band  structure  approach  it  is  natural  to  introduce  periodic  arrays 
of  non-interacting  defects  in  the  Si02  region  near  the  interface.  Before  embarking  on  such 
calculations,  we  first  performed  exploratory  studies  of  this  type  using  a  24-atom  Si02 
supercell  to  simulate  bulk  Si02,  as  reported  in  an  earlier  paper.*3  Similar  studies  were 
carried  out  some  time  ago  using  the  extended  Huckel  method.37  The  essential  idea  is  to 
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introduce  one  defect  such  as  a  Si  or  O  vacancy  or  interstitial  per  supercell,  and  then  calcu¬ 
late  the  energy  band  structure  at  two  or  more  well-separated  points  in  the  reduced  zone.  In 
this  way  one  obtains  not  only  the  valence  and  conduction  bands  of  the  host  Si02,  but  also 
impurity  bands  corresponding  to  the  defects  introduced.  If  the  supercell  is  sufficiently  large, 
and  the  defects  sufficiently  localized,  the  impurity  bands  will  have  negligible  dispersion,  and 
the  centers  of  gravity  of  these  impurity  bands  will  describe  the  energy  levels  of  isolated 
defects. 

An  important  advantage  of  the  periodic  defect  array  approach  over  molecular  cluster 
methods  is  that  the  energy  levels  of  the  defects  are  located  relative  to  the  valence  and 
conduction  band  edges  of  the  host  crystal.  We  will  now  briefly  recapitulate  our  earlier  Si02 
defect  investigations.13 

By  trial  and  error  we  find  that  a  unit  cube  of  S102  containing  24  atoms  places  adjacent 
defects  sufficiently  far  apart  (7.68  A)  to  produce  impurity  bands  with  dispersions  of  the 
order  of  0.1  eV,  which  is  small  enough  for  our  purposes.  For  some  defects  it  is  even  possible 
to  use  a  smaller  unit  cell  (Si02  with  6  atoms  per  unit  cell)  and  still  end  up  with  acceptably 
small  impurity  band  dispersions.  To  illustrate  our  approach,  we  consider  the  24  atom 
supercell,  which  is  constructed  by  placing  Si  atoms  at  (0,0,0),  (1/4, 1/4, 1/4),  etc.,  and  O 
atoms  at  (1/8, 1/8, 1/8),  etc.,  where  distances  are  measured  in  units  of  the  cube  edge.  The 
intrinsic  Si02  supercell  contains  128  valence  electrons,  so  there  are  normally  64  filled  bands. 
We  test  the  dispersion  properties  of  defect-related  bands  by  determining  the  band  structure 
of  the  supercell  at  the  zone  center  and  also  at  the  zone  corner.  The  essential  results  are  as 
follows: 

5.1  Interstitial  Silicon  Atoms 

Let  us  introduce  a  neutral  silicon  atom  at  the  interstitial  position  located  at  the  center  of  the 
Si02  cage  in  the  Si02  supercell,  (1/2, 1/2,  1/2).  The  64  filled  valence  bands  are  perturbed 
only  slightly.  There  is  a  filled  interstitial  band  roughly  midway  in  the  Si02  band  gap  arising 
from  the  Si  3s  level,  containing  2  of  the  4  valence  electrons  contributed  to  the  Si02 
supercell  by  the  Si  interstitial.  There  are  also  empty  interstitial  bands  located  well  above  the 
Si02  conduction  band  minimum  (CBM).  The  lowest  (perturbed)  conduction  band  is  also 
filled,  accounting  for  the  remaining  2  electrons  associated  with  the  Si  interstitial. 

Taking  into  account  the  nature  of  a  band-theoretic  description  of  a  periodic  array  of 
non-interacting  defects  in  an  otherwise  perfect  crystal,  the  above  results  can  be  interpreted 
as  follows:  If  we  introduce  a  Si +4  Interstitial  ion  core  into  Si02,  plus  the  4  electrons 
required  to  establish  charge  neutrality,  2  of  these  electrons  will  occupy  a  level  in  the 
mid-range  of  the  forbidden  band,  while  the  other  2  will  be  ionized,  occupying  low-lying 
conduction  band  states.  This  view  is  somewhat  naive,  reflecting  the  limitations  of  a  model  in 
which  the  defects  are  separated  just  far  enough  so  that  they  don’t  interact  to  any  appreciable 
degree.  In  practice,  the  final  2  electrons  required  for  local  charge  neutrality  would  not  be 
uniformly  distributed  throughout  the  crystal,  but  would  be  concentrated  on  the  atoms 
defining  the  cage  surrounding  the  Si  interstitial  atom. 

It  is  interesting  to  find  that  a  Si  interstitial  produces  a  level  roughly  midway  in  the  Si02 
forbidden  band.  This  suggests  that  Si  interstitials,  if  they  actually  exist  in  Si02,  could  create 
levels  within  striking  distance  of  the  thermal  gap  of  the  Si/Si02  interface.  Our  calculations 


Appendix  B 


Silicon-Silicon  Dioxide  Interfaces 


36 


are  too  crude  at  this  stage  and  the  possible  positions  of  the  Si  interstitial  too  numerous  in 
actual  Si02  for  us  to  make  predictions  more  precise  than  this  just  now. 

In  a  subsequent  study,14  we  investigated  the  effect  of  putting  a  neutral  interstitial  Si  atom  at 
the  center  of  a  Si02  cage  lying  next  to  a  Si/Si02  interface.  In  this  particular  study  we 
added  the  interstitial  Si  atom  to  a  17  atom  Si/Si02  superlattice  (cf.  Fig.  1  in  Ref.  12) 
composed  of  5  Si  layers  (2  atoms  each),  one  tier  of  Si04  tetrahedra  (5  atoms),  and  2  O 
atoms  doubly  bonded  to  the  two-fold  coordinated  Si  atoms  at  the  interfaces.  For  this 
particular  geometry,  the  interstitial  Si  atom  is  close  enough  to  the  doubly  bonded  O  atoms  to 
interact  with  them,  leading  to  a  lowering  of  the  Si  3s  and  3p  atomic  levels,  relative  to  their 
positions  in  the  Si02  super  cell.  The  net  result  is  that  the  Si  3s  band  lies  well  below  the 
Si/Si02  thermal  gap,  while  the  Si  3p  bands  lie  within  this  gap. 

These  studies  renfesent  the  first  stage  of  an  investigation  of  interstitial  Si  atoms  at  an 
interface  interacting  with  O  atoms  originally  doubly  bonded  to  the  silicon  substrate.  The 
object  is  to  study  the  formation  of  Si-O-Si-O-Si  bridges  in  which  the  O  atoms  are  singly 
bonded  to  the  Si  substrate  atoms  as  well  as  to  the  bridging  Si  atom  (originally  the  intersti¬ 
tial).  The  bridging  Si  atom  could  be  doubly  charged,  representing  one  form  of  fixed  positive 
charge  at  the  interface. 


5.2  Interstitial  Oxygen  Atoms 

We  have  also  introduced  neutral  O  interstitial  atoms  at  the  center  of  one  of  the  cages  in  the 
24-atom  S102  supercell,  position  (1/2, 1/2, 1/2).  This  leads  to  an  occupied  interstitial  O  2s 
band  about  1 1  eV  below  the  Si02  valence  band  maximum  (VBM),  two  occupied  O  2p  bands 
about  2  eV  above  the  VBM,  and  an  empty  O  2p  band  about  3  eV  above  the  VBM.  In 
practice,  such  interstitial  O  atoms  would  be  expected  to  complete  their  valency  by  attaching 
themselves  to  a  pair  of  H  atoms,  for  example,  to  form  water. 

If  we  place  the  neutral  O  interstitial  atom  at  the  center  of  the  Si02  cage  in  the  17-atom 
Si/Si02  superlattice  described  above,  we  find  that  the  O  2s  and  2p  impurity  bands  lie  about  7 
to  8  eV  below  their  corresponding  values  in  the  Si02  supercell.  The  reason  again  is  the 
proximity  of  the  interstitial  O  to  the  nearby  O  atoms  attached  to  the  silicon  substrate.  This 
situation  represents  a  step  in  the  direction  of  creating  a  Si-O-O-O-Si  bridge,  where  the 
interstitial  O  lies  at  the  middle  of  the  span.  Because  the  O  interstitial  has  acquired  its  full 
valency  of  8  electrons,  the  highest  valence  band  is  empty,  simulating  a  shallow  acceptor  level. 

These  studies  suggest  that  the  O  interstitial  is  not  lively  to  produce  localized  levels  within 
the  range  of  the  thermal  (Si)  gap,  in  contrast  to  the  Si  interstitial.  Roughly  speaking,  the 
highest  occupied  interstitial  O  levels  will  tend  to  line  up  with  the  highest  occupied  O  levels  in 
Si02  (within  a  few  eV),  this  energy  range  being  well  below  the  thermal  gap. 

5.3  Neutral  Si  Vacancies 

Returning  to  the  normal  24-atom  Si02  supercell,  let  us  now  introduce  a  neutral  Si  vacancy 
by  removing  a  Si  atom,  say  from  position  (0,0,0).  This  removes  4  electrons  from  the 
supercel),  so  the  lowest  62  rather  than  the  lowest  64  valence  bands  are  occupied.  Because 
the  uppermost  valence  bands  in  the  normal  supercell  are  built  up  from  O  2p  non-bonding 
orbitals,  the  removal  of  a  Si  atom  has  negligible  influence  on  the  energies  of  these  bands.  As 
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a  consequence,  the  2  vacant  acceptor  levels  (bands  63  and  64)  remain  essentially  degenerate 
with  the  topmost  valence  band  (band  62).  If  these  acceptor  levels  become  occupied,  a 
negatively  charged  Si  vacancy  would  be  created.  The  situation  remains  essentially  the  same 
if  we  remove  a  Si  atom  from  the  center  of  the  Si04  tetrahedron  in  the  17-atom  Si/Si02 
superlattice. 

However,  if  in  this  superlattice  we  remove  the  four-fold  coordinated  Si  atom  at  the  inter¬ 
face,  i.e.,  the  atom  common  to  the  Si  and  Si02  regions,  the  dominant  effect  is  the  creation 
of  two  dangling  Si  bonds  and  the  concomitant  formation  of  localized  electronic  states  at  the 
interface.  The  levels  associated  with  the  non-bridging  O  atoms  are  similar  to  those  generat¬ 
ed  by  the  removal  of  a  Si  atom  from  the  center  of  a  Si04  tetrahedron  far  away  from  the 
interface  (as  in  the  24-atom  Si02  supercell). 

5.4  Neutral  Oxygen  Vacancies 

Returning  once  more  to  the  Si02  supercell  and  removing  a  neutral  O  atom,  we  are  left  with 
128  -  6  s  122  electrons  in  the  unit  cell.  These  electrons  fill  the  lowest  60  valence  bands  as 
well  as  one  O  vacancy  band  located  about  2  eV  below  the  CBM.  On  the  other  hand,  if  we 
remove  an  O  atom  from  the  Si-O-Si  bond  which  connects  directly  to  the  Si/Si02  interface, 
we  create  dangling  Si  bonds  and  localized  interface  states  in  the  thermal  gap.  In  all  the 
studies  reported  here,  we  have  not  attempted  to  relax  the  lattice  after  the  introduction  of  a 
vacancy  or  interstitial.  Some  of  the  most  interesting  defects  of  all  involve  vacancies  and 
interstitials  associated  with  relaxed  lattices9- 1 1  and  vacancies  leading  to  broken  bonds  which 
are  subsequently  saturated  by  the  attachment  of  H  or  OH  groups.9"11 


6.  INTERSTITIAL  02,  SiO,  AND  H20  MOLECULES 

We  have  also  placed  neutral  O7,  SiO,  and  H20  molecules  inside  the  Si02  cages  in  the 
17-atom  Si/Si02  superlattice.1'*  Each  molecule  was  inserted  so  that  its  center  of  gravity 
coincides  with  the  cage  center,  and  its  longest  dimension  is  perpendicular  to  the  interface. 
For  02,  which  normally  has  12  valence  electrons,  we  find  that  there  are  7  occupied 
02-related  bands  (2  electrons  of  opposite  spin  per  band).  Our  calculations  thus  indicate  that 
the  02  molecule  has  acquired  2  additional  electrons,  becoming  doubly  negatively  charged. 
The  highest  valence  band  is  now  empty,  corresponding  to  a  shallow  acceptor  level.  We  are 
presently  investigating  the  possibility  that  02  molecules  close  to  an  interface  can  form 
bridges  between  adjacent  pairs  of  two-fold  coordinated  Si  substrate  atoms.  To  form 
Si-0»0-Si  bridges,  the  terminal  Si  atoms  would  have  to  contribute  one  dangling  bond  each 
to  the  spanning  0«0  molecule. 


For  SiO,  which  normally  has  10  valence  electrons,  we  find  that  there  are  six  occupied 
SiO-related  bands,  so  the  SiO  molecule  acquires  2  additional  electrons,  emptying  the  highest 
valence  band  and  producing  a  shallow  acceptor  level.  For  H20,  which  normally  has  8  valence 
electrons,  there  are  4  occupied  H20-re!ated  bands,  so  the  molecule  remains  neutral.  Of 
course,  these  results  apply  to  the  particular  interstitial  position  chosen,  but  the  conclusions 
are  instructive  nevertheless. 


Appendix  B 


Silicon-Silicon  Dioxide  Interfaces 


38 


7.  SUBSTITUTIONAL  A1  AND  SUBSTITUTIONAL  A1  +  INTERSTITIAL  H 

Finally,  let  us  consider  the  substitution  of  a  neutral  A1  atom  for  a  neutral  Si  atom  in  the 
S102  supercell.  The  valence  and  conduction  bands  of  S102  are  not  affected  significantly  by 
this  replacement.  The  removal  of  one  electron  from  the  supercell  leads  to  a  half-empty 
highest  valence  band  (the  64  highest  band).  Within  the  present  framework,  this  represents  a 
partially  filled  acceptor  level  located  slightly  above  the  VBM  of  Si02.  We  can  fill  this 
acceptor  level  by  supplying  an  additional  electron  to  the  superceil.  For  purposes  of  illustra¬ 
tion,  let  us  introduce  a  neutral  interstitial  H  atom  close  to  the  substitutional  A1  atom  in  the 
Si02  supercell,  at  position  (-l/8,-l/8,-l/8).  The  H  atom  donates  its  electron  to  the  A1 
defect,  completing  the  occupation  of  all  the  valence  bands,  and  producing  an  empty  H  level 
roughly  midway  within  the  Si02  forbidden  band.  The  exact  position  of  the  empty  level 
depends  on  the  relative  positions  of  the  A1  and  H  atoms,  but  the  midrange  of  the  Si02  band 
gap  is  the  indicated  range.  Thus  we  see  that  such  a  complex  could  produce  localized  states 
within  the  thermal  gap  of  Si/Si02,  which  also  lies  in  the  midrange  of  the  SI02  band  gap. 
Similar  effects  would  be  expected  if  we  introduced  an  interstitial  monovalent  atom  such  as 
Na  instead  of  H. 


8.  SPECULATIONS  ON  THE  U-SHAPED  CONTINUUM  OF  INTERFACE  STATES 

Various  authors38  have  called  attention  to  the  fact  that  there  is  often  a  sharp  peak  in  the 
density  of  Si/Si02  interface  states  at  about  0.3  eV  above  the  Si  VBM,  and  that  this  peak 
can  be  removed  by  exposing  the  interface  to  hydrogen.  Once  this  peak  is  removed,  the 
distribution  of  interface  states  in  the  thermal  gap  of  Si  appears  to  be  U-shaped,  rising 
toward  the  valence  and  conduction  band  edges  and  falling  to  minimal  values  in  between.  A 
plausible  explanation38  for  the  disappearance  of  the  sharp  peak  is  that  the  dangling  Si  bonds 
at  the  interface  responsible  for  this  peak  are  removed  by  the  formation  of  Si-H  bonds. 

The  origin  of  the  U-shaped  continuum  is  not  understood.  There  are  obviously  many  possible 
explanations.  Using  our  own  studies  as  a  guide,  and  in  particular  the  fact  that  substitutional 
Ai  +  interstitial  M  can  lead  to  interface  states  in  or  near  the  thermal  gap  (M  *»  monovalent 
metal  atom  or  H),  we  will  go  one  step  further  and  ask:  Are  there  different  Si  positions  at  or 
near  the  interface  at  which  substitutional  Al  +  interstitial  M  could  produce  different 
interface  state  distributions,  for  example,  one  peaked  below  the  Si  VBM  and  another  above 
the  Si  CBM,  so  that  the  superposition  of  the  tails  from  these  two  peaks  would  resemble  a 
U-shaped  curve  in  the  thermal  gap? 

The  explanation  we  propose  is  that  there  are  Al  impurities  near  the  interface,  some  occupy¬ 
ing  the  centers  of  A104  tetrahedra,  and  others  the  four-fold  coordinated  positions  at  the 
interface  where  the  Si  and  S102  regions  are  joined.  The  latter  Al  impurities  occupy  the 
centers  of  distorted  Si2-Al-02  tetrahedra.  Allowing  for  the  different  positions  that  neigh¬ 
boring  interstitial  M  atoms  would  occupy  at  an  actual  interface,  we  would  get  two  distinct 
distributions  of  interface  levels,  one  offset  in  energy  from  the  other.  The  distribution 
corresponding  to  A104  +■  M  would  lie  above  the  one  corresponding  to  Si2-A!-02  +  M,  so 
the  tail  from  the  former  would  account  for  the  upper  part  of  the  U,  and  the  tail  from  the 
latter  for  the  lower  part.  The  individual  pieces  of  this  explanation  are  consistent  with  our 
calculations.  Finally,  we  note  that  the  exposure  of  the  interface  to  H  can  be  regarded  not 
only  as  a  means  for  saturating  the  dangling  Si  bonds,  but  also  as  a  means  for  forming  the 


Appendix  B 


Silicon-Silicon  Dioxide  Interfaces 


39 


AIO4  +  H  and  Si2-AI-C>2  +  ^  complexes  that  could  be  responsible  (in  whole  or  in  part)  for 
the  U-shaped  continuum. 


9.  CONCLUDING  REMARKS 

We  have  shown  how  idealized  Si/Sl02  superlattices  can  be  used  to  describe  the  Si/S102 
interface  and  to  construct  simple  models  of  defects  such  as  SI  and  O  vacancies  and  intersti¬ 
tials.  The  calculations  we  have  discussed  represent  early  attempts  at  dealing  with  complex 
geometrical  situations.  The  results  to  date  have  been  sufficiently  encouraging  to  suggest 
that  a  more  detailed  understanding  of  these  and  other  defects  can  be  achieved  by  carrying 
out  first  principles  calculations  self -consistently,  taking  lattice  relaxation  effects  into 
account. 

The  author  is  grateful  to  many  individuals  for  stimulating  discussions  and  correspondence, 
particularly  R.  S.  Bauer,  B.  E.  Deal,  D.  J.  Henderson,  B.  E.  Hobbs,  N.  M.  Johnson,  R.  V. 
Kasowski,  A.  G.  Revesz,  and  W.  E.  Rudge. 
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FIGURE  1.  Si/Si02  superlattices  and  interfaces.  Panel  (a)  shows  a  view  of  the  upper 
interface  and  the  Si02  region  as  seen  from  just  below  the  Si02  region  in  panel  (b).  Panels 
(b)  and  (c)  represent  front  views  of  superlattices  having  17  and  23  atoms  per  unit  cell, 
respectively.  Both  superlattices  have  five  silicon  layers  (2  atoms  each).  All  the  two-fold 
coordinated  Si  interface  atoms  are  doubly  bonded  to  O  atoms.  In  panel  (b)  there  is  one  tier 
of  Si04  tetrahedra  (5  atoms),  while  in  panel  (c)  there  are  three  such  tiers  (11  atoms).  The 
numbers  denote  the  planes  in  which  the  various  atoms  lie,  the  2  representing  the  front-most 
plane,  etc.  For  clarity,  the  O  atoms  in  the  backmost  (-2)  plane  are  not  shown.  As  best  seen 
in  panel  (a),  the  interface  Si  atoms  (heavier  outer  circles)  form  a  checkerboard,  with  the  Si 
atoms  attached  to  the  Si04  tetrahedra  corresponding  to  "red"  squares,  and  the  Si  atoms 
doubly  bonded  to  O  atoms  corresponding  to  "black"  squares.  The  "bulls-eyes"  in  panel  (a) 
represent  the  Si»0  groups  as  seen  from  the  O  end.  For  another  view  of  these  superlattices, 
see  C.  R.  Helms,  J.  Vac.  Sci.  Technol.  Id,  608  (1979). 
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ABSTRACT:  We  discuss  some  new  ideas  for  improving  first-principles 

muffin-tin-orbital-type  electronic  structure  calculations  for  complex  crystals  such  as  ternary 
and  multinary  compounds.  The  new  ideas  include:  construction  of  generalized  muffin-tin 
orbitals  and  energy-dependent  double  MTO  basis  sets;  cellular  decomposition  based  on 
strongly  overlapping  atomic  spheres;  and  evaluation  of  matrix  elements  in  the  residual 
interstitial  region  by  extrapolation. 


In  the  mid-1970s,  Andersen  and  his  coworkers1  developed  a  simplified  version  of  the  linear 
combination  of  muffin-tin  orbitals  (LCMTO)  method  based  on  the  atomic  sphere 
approximation  (ASA).  This  simplified  approach  has  been  used  successfully  to  studv  the 
electronic  structure  of  many  types  of  close-packed  solids,  particularly  transition  metals.*  The 
original  LCMTO  method3  has  by  now  been  superseded  by  the  more  sophisticated  extended 
muffin-tin  orbital  (EMTO)  method.4  In  a  recent  paper,5  we  discussed  extensions  of  the 
LCMTO/EMTO  method  aimed  at  increased  accuracy  and  reduced  running  times.  Our  own 
work  focuses  on  loosely  packed  structures  containing  many  atoms  per  unit  cell,  including 
surfaces  and  interfaces  represented  by  repeating  slab  geometries.5,5  In  the  present  paper  we 
will  indicate  the  essential  ideas  underlying  these  extensions.  These  ideas  are  well  suited  to 
first-principles  studies  of  ternary  and  multinary  compounds  such  as  chalcopyrites  and  their 
derivatives. 

It  is  possible  to  apply  the  ASA  to  some  types  of  loosely  packed  solids  after  first  transforming 
them  into  quasi-close-packed  solids  by  replacing  part  of  the  interstitial  region  with  "empty" 
spheres.  In  this  way,  most  of  the  physically  relevant  space  is  spanned  by  weakly  overlapping 
atomic  and  "empty"  spheres.  However,  this  transformation  is  practical  only  for  certain 
geometries,7  and  even  for  these  geometries  the  adequacy  of  this  approach  can  be  obscured  by 
the  arbitrary  choice  of  "empty"  sphere  radii. 

For  more  general  types  of  open  structures,  it  becomes  necessary  to  use  more  rigorous 
treatments  such  as  LCMTO/EMTO  methods,  which  unfortunately  require  significantly 
longer  running  times.  In  these  methods,  the  wave  functions  are  represented  by 
exponentially-damped  MTOs  (expMTOs)  with  tails  having  the  form  exp(-xr)/r,  as  well  as  by 
oscillatory-damped  MTOs  (oscMTOs)  with  tails  having  the  form  cos(«cr)/r.  Evaluation  of 
matrix  elements  in  the  interstitial  region  by  three-dimensional  integration  accounts  for  the 
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greatly  increased  computational  effort 

Although  first-principles  LCMTO/EMTO  calculations  hare  been  carried  out  successfully  for 
many  different  types  of  surfaces.  Interfaces,  and  loosely  packed  solids,  some  having  as  many 
as  44  atoms  per  unit  cell,®’®  the  usefulness  of  these  calculations  could  be  greatly  enhanced  if 
their  accuracy  could  be  increased  without  substantially  increasing  the  running  time.  The 
same  would  be  true  if  the  running  time  could  be  decreased  by  one  or  two  orders  of 
magnitude,  even  at  the  cost  of  a  slight  decrease  in  accuracy. 

In  order  to  improve  the  accuracy  of  LCMTO/EMTO  calculations,  it  is  necessary  to  increase 
the  variational  flexibility  of  the  MTO  basis  sets.  This  flexibility  is  presently  limited  by  the 
use  of  single-term  tail  functions,  such  as  exp(-xr)/r,  outside  the  MT-spheres;  and  by  the 
use  of  energy-  independent  radial  functions  inside  the  MT-spheres.  The  need  for  increased 
accuracy  is  particularly  acute  for  inhomogeneous  systems  such  as  surfaces  and  interfaces, 
where  localized  surface  and  interface  states  must  be  represented  by  compact,  highly 
adaptable  basis  functions.  The  single  expMTOs  presently  in  use  often  prove  inadequate  for 
certain  systems.  If  one  represents  these  states  by  a  linear  combination  of  two  expMTOs 
having  different  values  of  a,  as  suggested  by  Harris  and  Painter,®  the  variational  flexibility 
can  be  increased,  but  the  computational  effort  is  roughly  doubled.  The  reason  is  that 
modified  spherical  Bessel  functions  must  now  be  evaluated  for  two  values  of  a  rather  than 
just  one. 

The  computational  effort  is  also  roughly  doubled  if  a  double  MTO  basis  set  composed  of 
oscMTOs  as  well  as  expMTOs  is  introduced  in  an  attempt  to  improve  the  variational 
adaptability  of  the  basis  set.  This  approach  usually  works  well  for  homogeneous  systems  but 
often  fails  for  inhomogeneous  ones.  The  most  serious  difficulty,  at  least  for  some  types  of 
inhomogeneous  systems,  is  that  the  oscMTOs  are  of  longer-range  than  the  expMTOs, 
leading  to  undesirable  coupling  between  adjacent  surfaces  or  interfaces. 

We  have  enhanced  the  variational  flexibility  of  localized  MTO  basis  functions  by 
constructing  generalized  MTOs  composed  of  linear  combinations  of  single  expMTOs  and 
their  first  and  second  derivatives  with  respect  to  k.  Since  these  derivatives  can  be  expressed 
in  terms  of  the  leading  expMTO  by  recursion  relations,  and  since  the  generalized  MTOs 
satisfy  expansion  theorems  analogous  to  those  of  ordinary  expMTOs,  the  standard 
LCMTO/EMTO  formalism  is  readily  generalized.®  Although  the  generalized  formalism  is 
more  complicated  as  a  result  of  the  additional  variational  flexibility  built  into  the  basis 
functions,  the  computational  effort  is  expected  to  increase  only  slightly  over  that  required 
for  state-of-the-art  LCMTO/EMTO  calculations.  The  reason  is  that  ail  the  modified 
spherical  Bessel  functions  appearing  depend  only  on  a  single  value  of  «r,  and  all  these 
functions  can  be  generated  efficiently  with  the  aid  of  recursion  relations. 

Moreover,  energy-dependent  double-MTO  basis  sets  can  be  constructed  by  suitably 
rearranging  and  manipulating  the  matrix  element  information.  Since  this  procedure  will 
normally  lead  to  large  and  small  MTO  components,  it  is  possible  to  achieve  further 
computational  economies  by  solving  the  double-dimension  secular  equations  by  the  Lowdln 
partitioning  scheme.9  In  addition  to  improving  the  variational  adaptability  still  further,  the 
use  of  energy-dependent  basis  sets  enables  us  to  solve  the  secular  equations  over  a  wider 
energy  range  than  would  otherwise  be  possible. 
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In  order  to  make  the  study  of  extremely  large,  loosely  packed  systems  practical,  it  Is 
essential  to  reduce  the  computational  effort  drastically.  This  can  be  accomplished  only  by 
drastically  simplifying  the  most  time-consuming  step  in  LCMTO/EMTO  calculations, 
namely,  the  evaluation  of  the  interstitial  matrix  elements.  Instead  of  dealing  with  the 
interstitial  region  as  a  whole,  or  partitioning  this  region  into  space-filling  Wlgner-Seitz 
atomic  polyhedra,  as  in  the  original  LCMTO  method,5  we  surround  each  (non-overlapping) 
MT-sphere  with  a  larger,  concentric  outer  sphere,  and  then  replace  this  set  of  overlapping 
outer  atomic  spheres  by  an  equivalent  set  of  truncated  atomic  spheres  Limiting  the  radii  of 
the  outer  spheres  so  that  no  two  of  these  spheres  overlap  inside  a  third,  and  representing 
each  of  the  truncated  atomic  spheres  by  a  three-dimensional  Heaviside  operator,5,1®  the 
matrix  elements  of  this  operator  can  be  determined  analytically,  which  is  considerably  easier 
to  do  than  integrating  over  a  three-dimensional  mesh  as  would  be  necessary  for  polyhedral 
atomic  cells.  Moreover,  the  matrix  elements  themselves  can  be  determined  inside  the 
truncated  atomic  spheres  by  angular  and  radial  integrations,  which  is  again  considerably 
easier  than  three-dimensional  integration. 

In  earlier  scattered  wave  calculations  for  large  planar  organic  molecules,11  we  studied  the 
convergence  of  solutions  based  on  sequences  of  overlapping  atomic  sphere  models  having 
progressively  larger  radii.  Based  on  this  earlier  experience  with  open  structures,  we 
developed  an  extrapolation  procedure  for  obtaining  accurate  estimates  of  the  matrix 
elements  in  the  residual  interstitial  region,  Le.,  in  the  region  outside  the  overlapping  atomic 
spheres.5  The  essential  idea  is  to  use  overlapping  atomic  spheres  which  are  not  unreasonably 
large  for  the  angular  and  radial  integrations,  and  then  to  simulate  the  use  of  still  larger 
overlapping  atomic  spheres  by  extrapolation.  Of  course,  there  will  be  a  slight  loss  In 
accuracy  due  to  the  extrapolation,  but  this  should  be  insignificant  compared  with  the  gain  in 
accuracy  arising  from  the  use  of  Improved  MTO  basis  sets.  Although  there  is  some 
arbitrariness  in  the  choice  of  the  overlapping  sphere  radii,  the  consequences  of  this 
arbitrariness  can  be  virtually  eliminated  by  employing  carefully  designed  extrapolation 
procedures.5 

A  complete  description  of  the  mathematical  formalism  appears  in  Ref.  5,  and  a  detailed 
discussion  of  three-dimensional  Heaviside  operators  is  given  in  Ref.  10.  Detailed 
calculations  based  on  these  ideas  are  currently  in  progress,  and  will  be  reported  elsewhere. 

*  Supported  in  part  by  the  Office  of  Naval  Research,  Contract  Number 
N00014-79-C-0814. 
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FIGURE  1.  Truncated  atomic  sphere  model  for  TCNQ  molecule  based  on  strongly 
overlapping  atomic  spheres.  The  centers  of  the  C,  N,  and  H  spheres  are  indicated  by  symbols 
such  as  C5.  The  overlapping  spheres  are  represented  by  the  dashed  lines,  and  the  truncated 
spheres  by  the  heavy  lines.  This  model  satisfies  the  multiple  overlap  restriction;  that  is  to 
say,  different  overlap  regions  do  not  overlap  one  another.  The  function  of  the  extrapolation 
procedure  described  in  the  text  is  to  include  more  and  more  of  the  residual  interstitial  region 
inside  a  sequence  of  overlapping  atomic  sphere  models  having  progressively  larger  radii.  In 
the  limit,  these  spheres  become  the  Van  der  Waals  spheres,  whose  envelope  describes  the 
physically  relevant  portion  of  the  interstitial  region.  For  further  details,  see  Refs.  5,  10,  and 
11. 


Note.  The  published  paper  includes  a  mathematical  appendix  which  has  not  been  included 
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It  is  possible  to  account  for  measured  Pd/Si  Schottky  barrier  heights  (0.71-0.73  eV)  by  assuming  that 
the  interface  between  Si  and  PtLSi  includes  a  hexagonal  Si  transition  region  induced  by  the  presence  of 
Pd  impurities.  With  the  Fermi  level  pinned  slightly  above  the  Si  valence  bead  edge  by  hydridfaeed  Pd-Si 
states  at  the  interface,  as  calculations  indicate,  the  Schottky  barrier  height  is  determined  primarily  by 
the  band  gap  of  hexagonal  Si,  which  is  about  0.85  eV.  Since  Pd  atoms  ate  larger  than  Si  atoms, 
isolated  substitutions!  and  interstitial  Pd  atoms  are  not  readily  incorporated  into  the  Si  lattice.  We 
believe  that  the  Pd  impurities  form  triangular  dusters  at  Si  vacancies.  Since  thaae  cinema  are  more 
easily  accommodated  in  hexagonal  than  in  cubic  Si,  the  Si  stacking  sequence  changes  from  cable  to 
hexagonal,  giving  rise  to  the  transition  region. 


1.  INTRODUCTION 

In  spite  of  extensive  research  on  rectifying 
meut-eemiconductor  contacts,  our  understanding  of  the 
factors  determining  Fermi  level  pinning  and  Schottky 
barrier  height*  is  still  rudimentary. 1-3  In  order  to  gain 
further  insight  into  the  nature  of  transition  metal  - 
semiconductor  interfaces,  we  decided  to  study  the  Pd/Si 
system,  for  which  a  wealth  of  experimental  information  is 
available.4*3  In  contrast  to  earlier  theoretical  calculations 
of  idealized  Pd/Si  compounds4  and  Pd  monolayers 
chemisorbed  on  the  Si(  1 11 )  surface,7  our  studies  are 
concerned  explicitly  with  PtLSi/Si  interfaces.  Using  a 
repeating  slab  geometry  and  the  first-principles  extended 
muffin-tin  orbital  method,  we  investigated  the  electronic 
structure  of  various  PdjSi/Si  superlattices  as  well  as  the 
effects  of  introducing  Pd  impurities  into  the  Si  substrate. 
Because  of  space  limitations,  we  will  confine  our 
attention  here  to  one  particular  set  of  ideas  that  have 
emerged  from  our  work.  The  details  of  our  calculations 
win  be  reported  elsewhere.* 

2.  PdjSi/Si  INTERFACES  AND  SUPERLATTICES 

One  of  the  major  obstacles  to  understanding 
metal -semiconductor  interfaces  is  our  limited  knowledge 
of  the  atomic  structure  of  such  interfaces.  In  order  to 
avoid  many  of  the  problems  associated  with 
compoaitionally  graded  interfaces  tad  other 
manifestations  of  structural  disorder,  we  decided  to  study 
epitaxial  systems.  Our  intsrest  in  the  PdjSi/Si  system 
was  reinforced  by  the  fact  that  PtLSi  can  be  grown 
epitaxially  on  a  Si  substrate.9*0  To  a  first 
approximation,  tbs  PtLSi/ Si  interface  can  be  mgsrded  as 
t  PdjSi  crystal  in  perfect  registry  with  s  Si  crystal.  For 
computational  reasons,  it  is  convenient  to  study 
PtLSi/Si(lll)  supsriattices  whose  interfaces  are 
sufndsntly  far  apart  to  be  noninteracting.11  As  can  be 
seen  from  Figure  1,  each  Si(lil)  layer  contributes  3  Si 
atoms  to  a  superiettice  sail  cad,  while  the  Si-rich  and 
Pd-rich  PdjSi  Inyen  contribute  3  Pd  +  2  Si,  and  3  Pd  + 
1  Si.  respectively. 


Figure  1.  The  Si(lll)  plane  is  shown  at  the  top, 
the  Si-rich  layer  of  PdjSi  in  the  middle,  and  the 
Pd-rich  layer  of  PdjSi  at  the  bottom.  In  the  top 
panel,  the  symbols  C,  B,  and  A  denote  the 
positions  of  Si  atoam  in  different  SKI II)  layers. 
We  also  use  these  symbols  to  denote  individual  Si 
layen  (not  double  SI  layers)  in  superiettice*. 
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83g  F.  Herman  etal.  /  Electronic  states  and  Schottky  barriers  at  Pd^SUSl  (111)  interfaces 


We  investigated  the  electronic  structure  of  various 
PdjSi/SKUl)  superiattices.  The  repeat  unit  of  one  of 
the  simplest  (8  atomic  layers)  contains  6  layers  of  cubic 
Si  (stacking  sequence  A-C-C-B-B-A)  followed  by  the 
Si-rich  and  Pd-rich  layers  of  PdjSi.  Altogether,  the  unit 
cell  contains  6  Pd  +  21  Si  atoms.  We  also  investigated 
an  even  simpler  6-layer  PdjSi/SKUl)  superlattice 
compoeed  of  4  layers  of  hexagonal  Si  (stacking  sequence 
(A-B-B-A)  and  the  Si-rich  and  Pd-rich  layers  of  Pt^Si. 
The  unit  ceil  here  contains  6  Pd  +  15  Si  atoms. 

For  all  such  superlattices,  there  are  clearly  two  different 
types  of  PdjSi/Si  interfaces,  corresponding  to  Si-rich  or 
Pd-rich  frontier  layers  of  PtLSL  For  both  types,  the 
relative  paeitkms  of  substrate  Si  atoms  and  adjacent  Pd 
atoms  are  different  from  the  on-top  and  triangular  sites 
«ttd  in  Ref.  7  to  describe  a  Pd  monolayer  chemisorbed 
on  the  Sl(  Ill)  surface.  Nevertheless,  using  reasonable 
estimates  of  Pd-Si  distances  across  the  interface  and 
more  realistic  interface  geometries,  our  own  preliminary 
results  are  consistent  with  those  of  Ref.  7  in  suggesting 
that  the  Fermi  level  is  pinned  by  hybridized  Pd-Si  states 
just  above  the  Si  valence  bend  edge.  The  pinning  is 
slightly  different  at  the  Si-rich  and  Pd-rich  PdjSi 
interfaces,  so  that  different  samples  having  a 
preponderance  of  one  interface  type  or  the  other  would 
have  shghriy  different  Schottky  barrier  heights.* 

3.  HEXAGONAL  SILICON  INDUCED  BY  Pd 
CLUSTERS 

In  order  to  account  for  the  fact  that  me  seined  PdjSi/Si 
Schottky  barrier  heights  (0.71-0.73  eV)5  are 
considerably  lower  than  theoretical  estimates  ( 1.0-1. 1 
eV),7  it  is  necessary  to  consider  deviations  from  ideally 
abrupt  interface  models  separating  perfect  PdjSI  and  Si 
crystals.  High  resolution  electron  microscope  studies10 
suggest  that  actual  Pdj Si/Si  interfaces  may  be  rough  on 
the  scale  of  a  several  atomic  layers,  may  consist  of  both 
interface  types,  and  may  contain  rotational ly  twinned 
domains  as  well.  Instead  of  exploring  such  complicated 
situations,  which  is  dearly  beyond  our  scope,  we  decided 
to  explore  the  effects  of  introducing  simpler  types  of 
defects,  including  Pd  impurities  and  stacking  disorder, 
into  our  Si  substrates.  Using  half  the  nearest  neighbor 
distances  in  Pd  metal  (2.73  A)  and  Si  (2.33  A)  as 
measures  of  atomic  size,  the  Pd  atom  is  17  percent  larger 
ihm  the  Si  atom,  and  to  would  not  be  readily 
accommodated  at  either  substitutional  or  interstitial  Si 
sites.  This  would  be  true  regardless  of  the  local  Si 
stacking  sequence. 

In  trying  to  imagine  how  Pd  atoms  would  arrange 
themselves  in  the  Si  subetrate  if  they  diffused  into  this 
substrate  in  advance  of  a  growing  interface,  we  realized 
that  triangular  Pd  clusters  such  as  thorn  belonging  to  the 
Si-rich  PdjSi  layer  could  be  accommodated  quite  well  at 
9  vacancies  in  a  hexagonal  Si  envirotunent,  and  less  well 
at  Si  vacancies  in  a  cubic  SI  environment.  The  hexagonal 
aequenca  A-B-B-A  leads  to  prismatic  hexagonal  cavities, 
while  the  cubic  sequence  A-B-B-C  leads  to 


adamantine-like  cavities.  Using  bell-and-stick  models,  it 
is  easily  seen  that  three  adjacent  interstitial  cavities 
surrounding  a  Si  vacancy  provide  more  room  for 
triangular  Pd  clusters  in  hexagonal  than  in  cubic  Si. 

We  were  originally  led  to  the  idea  of  introducing 
hexagonal  Si  regions  into  our  Pt^Si/Si  superlattices  in 
order  to  simplify  our  calculations  by  dealing  with  unit 
cells  containing  as  few  as  21  rather  than  27  atoms. 
However,  once  we  realized  that  triangular  Pd  dusters 
could  be  accommodated  more  easily  in  hexagonal  Si  than 
in  cubic  Si,  we  decided  to  pursue  this  idea  on  its  own 
merits  and  to  report  this  aspect  of  our  work  here. 
Accordingly,  let  us  set  aside  the  aaodel  of  an  ideal 
Pd^Si/Si  interface,  and  consider  instead  a  more  realistic 
interface  where  individual  Pd  atoms  or  small  Pd  dusters 
diffuse  into  the  Si  region  somewhat  In  advance  of  the 
growing  PdjSi/Si  interface. 

We  believe  that  individual  Pd  atoms  are  unlikely  to  enter 
the  Si  lattice  substitutionally  or  interstitially  because  of 
their  large  size,  but  instead  would  form  triangular  dusters 
which  would  replace  individual  Si  atoms.  During 
interface  formation,  we  would  expect  the  Si  stacking 
sequence  in  the  neighborhood  of  the  interface  to 
fluctuate.  In  order  to  accommodate  the  triangular  Pd 
clusters,  the  hexagonal  Si  structure  would  eventually 
form.  Even  though  the  total  energy  .of  a  hexagonal  SI 
lattice  is  higher  than  that  of  a  cubic  Si  lattice,  the  energy 
difference  is  quite  small,  and  is  easily  compensated  by 
the  reduction  in  elastic  energy  sseociated  with  the 
improved  accommodation  of  the  triangular  Pd  clusters  at 
Si  vacancies.  It  would  probably  take  only  a  small 
number  of  such  clusters  to  influence  the  stacking 
sequence  over  an  extended  region. 

In  short,  we  believe  that  a  hexagonal  Si  transition  region 
would  be  generated  just  in  front  of  the  growing  PdjSi/Si 
interface.  Bearing  in  mind  that  the  Si-rich  Pd^Si  layer  is 
an  ordered  triangular  array  of  triangular  Pd  clusters  and 
Si  atoms  (cf.  Figure  1),  a  Si  layer  containing  triangular 
Pd  clusters  at  random  isolated  Si  vacancy  sites  can  be 
regarded  as  an  incipient  Si-rich  PdjSi  layer.  It  is  difficult 
to  estimate  the  width  of  this  transition  region,  but  this 
width  would  probably  be  comparable  to  the  scale  of  the 
interface  roughness,  that  is  to  my,  several  atomic  layers. 
Even  if  the  hexagonal  Si  transition  region  consists  of 
only  several  atomic  layers,  its  electronic  structure  would 
be  characteristic  of  hexagonal  Si,  and  the  Schottky 
barrier  at  the  Pd2Si/Si  interface  would  be  determined  in 
large  measure  by  the  bandgap  of  hexagonal  Si,  rather 
than  by  that  of  the  cubic  Si  substrate. 

4.  BAND  STRUCTURE  OF  HEXAGONAL  SILICON 

It  is  well  known  that  direct  band  gap  semiconductors 
such  as  ZnS  have  nearly  the  same  band  gap  in  the  cubic 
and  hexagonal  structures.  It  is  perhaps  less  well  known 
that  the  polytypes  of  indirect  band  gap  materials  such  as 
SiC  have  different  band  gaps  because  the  conduction 
band  Minima  are  shifted  to  different  positions  in  reduced 
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zoom  having  different  symmetric!. 1 2  For  example,  in 
going  from  cubic  (3C)  SiC  to  hexagonal  (2H)  SiC,  the 
conduction  band  edge  shift!  from  the  X  point  to  the  K 
point  in  the  corresponding  reduced  zones,  and  the 
indirect  band  gap  changM  from  2.3  to  3.3  eV.12  Clearly, 
these  effects  can  be  quite  large. 
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ABSTRACT.  We  have  generated  an  atomic-scale  model  of  the  (001)  interface  between 
amorphous  and  crystalline  tetrahedrally-coordinated  semiconductors.  The  construction  is 
based  on  original  computer  algorithms  which  transform  random  and  ordered  close-packed 
lattices  into  corresponding  covalent  networks.  Calculations  of  the  electronic  density  of 
states  indicate  the  presence  of  gap  states  in  the  amorphous  region  and  near  the  interface. 
These  gap  states  arise  from  structural  imperfections,  primarily  dangling  bonds,  which  occur 
naturally  in  this  fully  computer-generated  model,  as  do  voids  and  lattice  distortions.  The 
spatial  extent  of  the  interface  is  a  few  atomic  layers. 


I.  INTRODUCTION 

Although  a  great  deal  of  theoretical  effort  has  been  devoted  to  the  study  of  idealized 
interfaces  connecting  lattice-matched  crystalline  semiconductors,1  interfaces  between 
crystalline  (c)  and  amorphous  (a)  semiconductors  have  received  only  scant  attention.2  To  onr 
knowledge,  the  only  existing  atomic-scale  models  of  such  interfaces  are  small  (abont  100 
atom)  hand-built  clusters.3-3  In  order  to  avoid  subjective  influences  in  the  construction  of 
such  models,  we  have  developed  a  fully  automatic  technique  and  applied  this  to  the  study  of 
the  c-Si/a-Si  interface.  (We  will  use  Si  for  purposes  of  discussion,  even  though  onr  struc¬ 
tural  results  apply  in  principle  equally  well  to  C,  Si,  Ge,  and  grey  Sn.)  As  a  first  step,  we 
generate  Lennard- Jones  crystal/melt  interfaces  by  Monte  Carlo  simulation.  This  leads  to 
dense  configurations  which  are  then  transformed  into  adjoining  ordered  and  disordered 
covalent  networks,  thereby  creating  an  atomic  model  of  the  interface  between  c-Si  and  a -Si. 
The  geometrical  transformation  between  close-packed  and  covalent  structures*  is  similar  to 
those  used  earlier  in  generating  bulk  amorphous  Si  clusters  from  random  hard-sphere 
configurations.7'*  Finally,  the  local  electronic  density  of  states  near  the  interface  is 
determined  by  tight-binding /recursion  methods. 


II.  AMORPHOUS/CRYSTAL  SEMICONDUCTOR  INTERFACES 

We  have  applied  the  transformation  from  dose-packed  to  covalent  networks  to  the 
crystal/melt  interface  of  the  Lennard-Jones  system9  governed  by  the  interatomic  potential 
V(r)  m  4«((o/r)12  -  (o/r)*j,  where  e  and  a  are  the  units  of  energy  and  length,  respectively. 
Onr  starting  configuration  is  a  fee  lattice  composed  of  27  (OOl)-planes,  each  containing  50 
atoms.  Our  choice  of  parameters  is  based  on  coexistence  data10  for  the  Lennard-Jones 
system  at  the  temperature  kT«1.35e:  In  the  ordered  regions  (planes  1  to  5  and  23  to  27), 
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the  density  was  set  equal  to  ps=1.053o-3,  and  in  the  disordered  region  to  =0.964a"^,  the 
density  of  the  liquid  phase.  The  corresponding  interplanar  spacings  are  0.780a  and  0.852o. 
The  400  atoms  in  planes  1  to  4  and  24  to  27  are  then  held  fixed  at  their  fee  positions,  while 
the  positions  of  the  950  atoms  in  the  remaining  planes  are  rearranged  according  to  Monte 
Carlo  techniques,  with  periodic  boundary  conditions  imposed  in  the  x  and  y  directions. 

Once  the  crystal/melt/crystal  Lennard- Jones  system  has  been  created,  we  suppress  the 
periodic  boundary  conditions  and  construct  the  Delauney  or  simplicial  graph.  This  graph 
subdivides  space  into  distinct  non-overlapping  tetrahedra  whose  vertices  are  the  atomic 
positions.11  The  simplicial  graph  for  our  system  contains  about  6700  tetrahedra.  After 
eliminating  ill-shaped  tetrahedra,12  specially  designed  algorithms  are  used  to  select  a 
simplicial  subgraph  which  only  includes  tetrahedra  not  sharing  any  common  edge  with  one 
another.  This  subgraph  is  not  unique  because  the  order  in  which  tetrahedra  are  considered 
for  inclusion  or  rejection  may  influence  the  selection  process.  The  simplicial  subgraph  is 
then  decorated  by  placing  Si  atoms  at  all  the  vertices,  as  well  as  at  the  centroids  of  all  the 
tetrahedra.  Keeping  all  the  vertex  atoms  fixed,  the  positions  of  the  centroid  atoms  are  then 
relaxed.  The  final  position  of  each  interior  atom  is  determined  by  minimizing  the  Keating 
elastic  energy13  of  the  five-atom  cluster  composed  of  the  interior  atom  and  the  four 
surrounding  vertex  atoms.  Finally,  bonds  are  established  between  the  vertex  atoms  of  the 
residual  tetrahedra  and  their  respective  interior  atom*.7-8  Since  the  vertex  atoms  of  the 
residual  subgraph  can  be  shared  by  fewer  than  four  tetrahedra,  dangling  bonds  and  voids  are 
normally  present  in  the  structure.  This  process  leads  to  residual  subgraphs  containing  about 
800  tetrahedra  and  to  covalent  networks  containing  about  2,100  atoms.  Since  the  method 
transforms  the  fee  lattice  into  the  diamond  structure,  atoms  located  in  planes  1  to  4  and  24 
to  27  of  the  initial  fee  configuration  are  automatically  arranged  in  the  diamond  structure. 
The  choice  of  geometry  clearly  insures  the  creation  of  the  (001)  c-Si/a-Si  interface. 

The  entire  process  just  described  was  repeated  five  times  in  order  to  generate  five  independ¬ 
ent  final  configurations.  The  mean  Interatomic  distance  in  the  amorphous  phase  was  found 
to  be  0.70 lo,  with  a  r.m.s.  deviation  of  Q.O680.  The  corresponding  values  for  the  bond- 
angle  distribution  were  110°  and  12°.  The  Si-Si  distance  in  the  crystalline  phase  is  0.676 cr, 
so  the  crystalline  phase  is  slightly  more  dense  than  the  amorphous  phase.  This  result  is  due 
partly  to  the  presence  of  voids  and  dangling  bonds  in  the  amorphous  region,  and  partly  to  the 
initial  choice  of  the  interplanar  spacings  in  the  ordered  and  disordered  regions.  (We  used 
slightly  different  interplanar  spacings  Initially  in  order  to  insure  melting  of  the  disordered 
region.)  In  Fig.  1  we  show  some  key  results  obtained  by  averaging  the  results  for  the  five 
final  configurations.  The  connectivity  (average  number  of  bonds  per  atom)  is  less  than  four 
in  the  amorphous  region  because  of  the  presence  of  dangling  bonds.  The  dangling  bond 
density  remains  roughly  constant  throughout  this  region,  except  near  the  interfaces,  where 
the  average  coordination  number  increases. 


III.  ELECTRONIC  DENSITY  OF  STATES  OF  THE  c-Si/a-Si  INTERFACE 

In  Fig.  2  we  show  local  densities  of  states  (LDOS)  that  were  calculated  using  a  first- 
neighbor  tight-binding  Hamiltonian  whose  two-center  integrals  were  obtained  by  fitting  a 
recent  LCAO  band  structure  of  c-Sl.14  The  moments  of  the  LDOS  were  computed  using 
the  recursion  method.  Averages  over  the  five  configurations  were  then  taken.  Since  only  a 
small  number  of  moments  was  included,  the  energy  resolution  of  the  LDOS  is  relatively  low. 
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Nevertheless,  the  existence  of  states  in  the  gap  is  clearly  seen  in  Fig.  2.  The  states  in  the 
gap  are  due  primarily  to  dangling  bonds,  and  to  a  lesser  degree  to  lattice  distortion.  The  gap 
mte  density  increases  as  one  moves  away  from  the  crystalline  region  across  the  interface 
into  the  amorphous  region,  the  range  of  variation  being  a  few  atomic  spacings. 


IV.  CONCLUDING  REMARKS 

In  summary,  we  have  succeeded  in  constructing  atomic-scale  models  of  c-Si/a-Si  interfaces 
fully  automatically.  These  computer-generated  interf;  ?es  are  a  few  atomic  spacings  wide 
and  are  undoubtedly  more  realistic  in  this  respect  than  atomically  abrupt  c-Si/a-Si  (and 
c-Si/a-Si02)  interfaces  that  have  been  constructed  by  hand  in  the  past.3'5  The  atomic 
models  could  be  improved  by  refining  the  construction  so  as  to  reduce  the  number  of  voids 
and  dangling  bonds.  Further  (but  probably  only  slight)  improvements  could  also  be  realized 
by  relaxing  all  the  atoms  in  the  system,  and  by  retaining  periodic  boundary  conditions  at  all 
stages  of  the  work.  The  present  approach  can  be  readily  extended  to  describe  interfaces 
between  c-Si  and  hydrogenated  a-Sl  or  a-Si02,  as  well  as  to  other  types  of  covalently  bonded 
networks. 
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Fig.  1.  Profiles  of  coordination  number  and  atomic  density  (in  units  of  atoms  per  a3) 
averaged  over  five  configurations.  Lateral  boundary  effects  are  avoided  by  ignoring  atoms 
lying  near  the  surface.  Fee  planes  in  the  initial  configuration  are  labeled  1  through  27. 
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Fig.  2.  Local  electronic  density  of  states  averaged  over  five  configurations.  Panels  A 
through  H  correspond  to  the  intervals  identified  in  Fig.  1.  Periodic  boundary  conditions 
were  applied  only  in  the  z  direction.  Only  atoms  lying  far  away  from  the  lateral  surfaces 
were  considered  so  as  to  avoid  boundary  effects.  The  zero  of  energy  corresponds  to  the  top 
of  the  valence  band  of  c-Si  OTjg.). 
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ELECTRONIC  AND  MAGNETIC  STRUCTURE  OF  ULTRATHIN 
COBALT-CHROMIUM  MULTILAYERS* 

F.  Herman  and  P.  Lambin,**  IBM  Research  Laboratory,  San  Jose,  California  95193 
O.  Jepsen,  Max  Planck  Institute  for  Solid  State  Research,  Stuttgart,  W.  Germany 

ABSTRACT:  In  order  to  study  exchange  coupling  and  spin  distributions  at  atomically  abrupt 
ferromagnetic/antiferromagnetic  interfaces,  we  construct  lattice-matched  Co/Cr  multilayer 
models.  Here  we  consider  strained-layer  superiattices  composed  of  alternating  regions  of 
ferromagnetic  bcc  Co  and  antiferromagnetic  bcc  Cr  with  repeat  periods  ranging  from  4  to  8 
atomic  layers.  For  computational  simplicity,  Cr  is  represented  by  a  theoretically-stabilized 
commensurate  antiferromagnetic  spin  arrangement.  The  multilayer  spin  distributions  are 
determined  by  carrying  out  first-principles  self-consistent  spin-polarized  LMTO/ASA 
electronic  structure  calculations.  We  find  that  the  magnetic  properties  of  the  Co/Cr 
superiattices  are  dominated  by  the  ferromagnetic  Co  layers,  though  the  antiferromagnetic 
character  of  bcc  Cr  is  still  evident.  For  multilayers  containing  1  Co  layer  and  3  to  7  Cr 
layers  per  repeat  period,  there  is  only  one  stable  spin  arrangement  corresponding  to  ferro¬ 
magnetic  coupling  across  the  Co/Cr  interfaces.  For  superiattices  containing  thicker  Co 
regions,  e.g.,  3  Co  layers  and  5  Cr  layers,  there  are  two  distinct  spin  arrangements  corre¬ 
sponding  to  ferromagnetic  as  well  as  antiferromagnetic  coupling  across  the  Co/Cr  interfaces. 
The  implications  of  these  results  for  more  complicated  Co/Cr  multilayers  are  discussed. 

PACS  numbers:  75.10.Lp,  75.70.  -  i,  75.50.Cc,  75.30.Cr 


INTRODUCTION 

In  this  paper  we  study  exchange  coupling  and  spin  distributions  in  crystallographicaily 
coherent  ferromagnetic/antiferromagnetic  multilayers  (superiattices)  composed  of  alternat¬ 
ing  ultrathin  slabs  of  Co  and  Cr.  Since  these  elements  have  nearly  the  same  atomic  volumes 
(within  about  2  percent),  it  is  reasonable  to  use  a  strained-layer  superlattice  as  the  theoreti¬ 
cal  model.  In  such  a  model  all  atoms  lie  on  a  common  lattice  having  suitably  averaged 
dimensions.  We  will  focus  on  Cr-rich  superiattices  having  repeat  periods  ranging  from  4  to 
8  atomic  layers,  with  [001]-oriented  interfaces.  Since  the  crystal  structure  of  strained-layer 
superiattices  should  be  dictated  by  the  majority  constituent,  we  will  assume  that  all  Co  and  Cr 
atoms  in  our  Cr-rich  structures  lie  on  a  common  bcc  lattice.  Although  we  regard  these 
hypothetical  multilayers  primarily  as  theoretical  models,  we  note  that  bcc  Co  overlayers  have 
recently  been  grown  on  bcc  Cr  substrates,1  suggesting  that  the  laboratory  synthesis  of  Co/Cr 
superiattices  also  bears  looking  into. 

We  will  determine  the  spin  distributions  in  bcc  Co/Cr  superiattices  by  carrying  out  self- 
consistent  spin-polarized  electronic  structure  calculations  using  frozen  cores,  the  first- 
principles  LMTO/ASA  method,2*3  and  the  spin-density  functional  approximation  of  von 
Barth  and  Hediu.4  We  recently  carried  out  similar  studies  for  ordered  fee  superiattices 
containing  ferromagnetic  NijFe  (Permalloy)  and  antiferromagnetic  MnFe.3  In  a  forthcom- 
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lag  publication, 6  we  will  discuss  the  present  work  more  fully,  as  well  as  lattice-matched 
superlattices  composed  of  ferromagnetic  hep  Co  and  antiferromagnetic  bcc  Cr. 


THEORETICAL  AND  COMPUTATIONAL  CONSIDERATIONS 

It  is  not  expedient  to  carry  out  first -principles  calculations  for  spin-density  ware  (SDW) 
antiferromagnetic  Cr  because  of  the  long  wavelength  (21  lattice  spacings).7  Moreover,  it  is 
likely  that  SDW-antiferromagnetism  would  be  suppressed  in  ultrathin  Cr  slabs.  We  turn 
instead  to  the  commensurate  antiferromagnetic  (C-AF)  model,8  according  to  which  bcc  Cr 
consists  of  two  interpenetrating  simple  cubic  lattices,  with  the  atoms  on  the  two  lattices 
having  opposite  spins.  Although  C-AF-Cr  provides  only  a  rough  approximation  of  SDW- 
AF-Cr,  we  can  justify  using  the  C-AF  model  for  three  reasons: 

First  and  foremost,  the  calculations  are  simplified  considerably.  Even  with  the  C-AF  model, 
8-layer  superlattice  studies  require  substantial  computational  effort.  Second,  the  C-AF  state 
of  bulk  Cr  can  be  stabilized  experimentally  by  the  addition  of  certain  transition  metal 
impurities,9  suggesting  that  ultrathin  antiferromagnetic  Cr  slabs  could  be  similarly  stabilized. 
Third,  the  impurity-stabilized  C-AF-state  of  Cr  is  not  likely  to  be  suppressed  by  the  1 
percent  or  so  reduction  in  lattice  dimensions  that  Cr  would  suffer  during  the  formation  of 
Co/Cr  strained-layer  superlattices.  In  contrast,  the  SDW-AF-state  of  bulk  Cr  is  destroyed 
by  even  smaller  reductions  in  lattice  dimensions.18 

In  order  to  gain  some  experience  with  Co/Cr  superlattices,  we  first  carried  out  some 
preliminary  calculations  for  bulk  Co  and  Cr,  and  then  examined  several  composite  systems, 
always  using  the  C-AF  model  for  Cr.  Contrary  to  earlier  studies,8  we  found  that  sustained 
self-consistent  iteration  of  C-AF-Cr  leads  to  a  paramagnetic  rather  than  an  antiferromag¬ 
netic  state.  Because  of  the  extremely  slow  convergence  of  antiferromagnetic  calculations, 
some  earlier  workers  had  concluded,  incorrectly,  that  their  solutions  had  converged  to  the 
C-AF  state.  We  are  presently  developing  improved  LMTO/ASA  computer  programs  which 
we  intend  to  use  to  re-examine  the  question  of  the  C-AF  state  of  bulk  Cr. 

As  an  outgrowth  of  these  studies,  we  developed  a  heuristic  method  for  theoretically  stabiliz¬ 
ing  the  C-AF  model  of  Cr.  Since  this  treatment8  appears  to  work  quite  well  for  composite 
systems  such  as  Co/Cr  multilayers,  as  well  as  for  bulk  Cr,  we  decided  to  base  the  present 
study  entirely  on  this  stabilized  C-AF  model,  which  we  regard  as  the  theoretical  equivalent  of 
experimental  impurity-stabilization.9 

Having  thus  set  our  course,  we  proceeded  to  study  the  convergence  properties  of  the  spin 
distributions  as  a  function  of  the  number  of  mesh  points  used  in  the  three-dimensional 
Integrations  over  the  reduced  zone  in  k-space.  For  materials  having  high  symmetry  and  few 
atoms  per  unit  cell,  such  as  (bulk)  hep  Co,  bcc  Co,  and  bcc  Cr,  It  is  necessary  to  use  many 
hundreds  of  mesh  points  in  order  to  obtain  say  1  percent  accuracy.8  Because  of  the  lower 
(tetragonal)  symmetry  of  our  superlattices,  and  the  larger  number  of  atoms  per  unit  cell, 
achieving  such  accuracy  would  require  enormous  computational  effort.  Accordingly,  we  had 
to  reach  a  compromise  between  accuracy  and  computer  time.  We  found  that  a  72-point 
mesh  was  adequate  for  our  purposes.  In  order  to  be  able  to  compare  the  magnitudes  of  Co 
and  Cr  spin  moments  In  the  superlattices  with  their  corresponding  bulk  values  —  at  the  same 


Appendix  F 


Cobalt-Chromium  Multilayers 


59 


level  of  convergence  —  we  also  studied  all-Co  and  all-Cr  4-,  6-,  and  8-layer  bcc  superlattic¬ 
es. 


STRAINED-LAYER  COBALT/CHROMIUM  SUPERLATTICES 
We  investigated  the  following  structures: 

•  4-layer  superlattices  —  Co(I)Cr(III),  Co(IV),  and  Cr(IV) 

•  6-layer  superlattices  —  Co(I)Cr(V),  Co(VI),  and  Cr(VI) 

•  8-layer  superlattices  —  Co(I)Cr(VII),  Co(III)Cr(V),  Co(VIII),  and  Cr(VIII), 

* 

where  roman  numerals  denote  the  number  of  successive  Co  or  Cr  atoms  in  the  repeat  period. 
For  example,  the  repeat  period  of  Co(I)Cr(III)  is  Co(l)Cr(2)Cr(3)Cr(4),  where  the  arable 
numerals  identify  the  atomic  site  or  layer.  We  started  off  our  calculations  for  these  Co/Cr 
superlattices  using  a  variety  of  assumed  spin  distributions,  and  then  iterated  to  self- 
consistency.  We  also  tried  different  starting  distributions  to  see  whether  we  could  generate 
alternate  self-consistent  solutions.  The  only  restraint  we  imposed  was  that  all  atoms  in  the 
same  atomic  layer  have  the  same  spin  moment  and  orientation,  as  would  be  expected  for 
[0011-oriented  interfaces  and  the  C-AF  model  for  Cr.  The  magnetic  properties  of  all  the 
composite  superlattices  were  dominated  by  the  ferromagnetic  Co  layers. 

For  the  Co/Cr  superlattices  containing  only  one  Co  layer,  namely,  Cod)Cr(III),  Co(I)Cr(V), 
and  Co(I)Cr(VII),  there  is  only  one  self-consistent  spin  distribution,  and  this  has  the  follow¬ 
ing  key  features: 

•  The  spin  moments  of  the  interfacial  Cr  atoms  are  aligned  parallel  to  the  spin  moments 
of  the  Co  atoms  (ferromagnetic  coupling  across  the  interfaces). 

•  The  spin  moments  of  the  Cr  atoms  alternate  from  layer  to  layer  (antiferromagnetic 
arrangement). 

•  The  spin  moments  of  the  Co  atoms  are  reduced  and  the  spin  moments  of  the  interfacial 
Cr  atoms  are  enhanced  relative  to  their  bulk  values. 

Generalizing  these  results,  we  anticipate  that  the  interfaciai  coupling  will  always  be  ferro- 
i  magnetic  for  Co/Cr  superiattices  having  single  Co  layers  separated  by  multiple  Cr  layers.  In 

the  simplest  cases,  successive  Co  layers  will  have  parallel  or  antiparallel  spins  according  as  the 
number  of  intervening  Cr  layers  is  odd  or  even.  More  complicated  spatially  modulated  spin 
distributions  can  also  be  visualized. 

For  Co(III)Cr(V),  which  contains  3  adjacent  Co  layers  separated  by  5  adjacent  Cr  layers,  the 
essential  results  are  as  follows: 

•  There  are  two  distinct  spin  distributions,  corresponding  to  parallel  and  antiparallel 
alignments  of  the  interfaciai  Co  and  Cr  spins  (ferromagnetic  as  well  as  antiferromag¬ 
netic  interfaciai  coupling). 

•  For  both  cases,  the  Cr  regions  are  antiferromagnetic  (C-AF  model). 

•  For  both  cases,  all  Co  and  Cr  spin  moments  are  reduced  relative  to  their  respective  bulk 
values;  moreover,  the  Co  moments  for  the  interfaciai  layers  are  smaller  than  the 
moments  of  the  central  Co  layer. 
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Since  there  was  only  one  spin  distribution  for  Co(I)Cr(V),  which  also  contains  5  adjacent  Cr 
layers,  we  attribute  the  added  degree  of  freedom  in  Co(III)Cr(V)  to  the  presence  of  more 
than  1  adjacent  Co  layer.  The  added  flexibility  arises  from  the  possibility  of  moment 
redistribution  in  the  interfacia]  and  central  Co  layers.  We  would  expect  analogous  results  for 
still  thicker  Co  slabs. 

The  general  features  just  found  for  (001]  Co/Cr  superlattices  should  also  apply  to  [111] 
Co/Cr  superlattices,  because  all  atoms  in  interfacial  Cr  layers  have  the  same  spin  for  both  of 
these  orientations  (C-AF  model).  On  the  other  hand,  for  (110]  Co/Cr  superlattices,  the 
interfacial  Cr  layers  would  be  compensation  planes  (equal  numbers  of  up  and  down  spins). 
Since  our  calculations  indicate  that  only  ferromagnetic  coupling  is  supported  across  single¬ 
layer  Co  interfaces,  we  would  expect  the  C-AF  state  of  Cr  to  be  suppressed  near  [110] 
single-layer  Co  interfaces.  By  the  same  token,  we  would  expect  the  C-AF  state  as  well  as 
more  complicated  antiferromagnetic  spin  orderings  of  Cr  to  be  sustained  near  [110] 
multiple-layer  Co  interfaces. 


CONCLUDING  REMARKS 

To  recapitulate,  we  have  calculated  the  exchange  coupling  and  spin  distributions  in 
ferromagnetic/antiferromagnetic  strained-layer  superlattices  composed  of  ultrathin  bcc  Co 
and  bcc  Cr  slabs.  By  calling  attention  to  these  hypothetical  superlattices,  which  are  used 
here  as  convenient  theoretical  models,  we  hope  to  stimulate  interest  in  the  laboratory 
synthesis  of  such  artificially  layered  structures,  as  well  as  the  experimental  study  of  their 
electronic  and  magnetic  properties.  We  are  currently  exploring  the  possibility  of  incorporat¬ 
ing  magnetocrystalline  anisotropy  Into  the  present  theoretical  model, 6  with  a  view  to  carrying 
out  first-principles  investigations  of  exchange  anisotropy.1*  For  further  details  concerning 
the  present  study,  see  Ref.  6. 
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TABLE  1.  Spin  magnetic  moments  for  4-laver  bcc  superlattices  Co(IV),  Cr(IV),  and 
Co(I)Cr(III)  *  Co(l)Cr(2)Cr(3)Cr(4),  for  various  meshes.  For  Co(I)Cr(III),  atomic  sites 
Cr(2)  and  Cr(4)  are  equivalent.  Here  and  in  the  following  tables,  Cr  is  represented  by 
theoretically  stabilized  commensurate  antiferromagnetic  model  discussed  in  Ref.  6. 


Mesh  for  tetragonal 

Co(IV) 

CrdV) 

< - 

-  Cod)Cr(III)  — 

- > 

reduced  zone 

(ail-Co) 

(all-Cr) 

Co(l) 

Cr(2)  Cr(3) 

Cr(4) 

IK- 20  (72)(a> 

1.56 

±  0.43 

+  1.37 

+  0.88  -  0.44 

+  0.88 

(a)  IK  is  number  of  inequivalent  mesh  points  in  irreducible  sector  of  reduced  zone.  Total 
number  of  mesh  points  in  reduced  zone  is  shown  in  parentheses.  See  Ref.  6  for  a 
discussion  of  the  convergence  of  the  spin  moments  as  a  function  of  the  number  of  points 
in  the  mesh. 


TABLE  II.  Spin  magnetic  moments  for  bcc  superlattices  Co(VI),  Cr(VI),  and  Co(I)Cr(V) 
—  Co(l)Cr(2)Cr(3)Cr(4)Cr(5)Cr(6).  Atomic  sites  in  same  column  are  equivalent. 


Mesh  for  tetragonal 
reduced  zone 

Co(VI) 

(all  sites 
cobalt) 

Cr(VI) 

(all  sites 
chromium) 

< - 

Co(l) 

- Co(I)Cr(V)  - 

Cr(2)  Cr(3) 

CK6)  Cr(5) 

Cr(4) 

IK-20  (72) 

1.65 

±  0.38 

+  1.46 

+  0.94 

-0.42 

+0.36 

1 

I 

5  TABLE  III.  Spin  magnetic  moments  for  8-layer  bcc  superlattices  Cr(VIII)  and 

|  Cod)Cr(VII)  —  Co(l)Cr(2)  ...  Cr(8).  Atomic  sites  in  same  column  are  equivalent. 

Cr(VIII) 
(ail  sites 
chromium) 

Co(I)CKVH) 

Cr(3) 

Cr(7) 

Mesh  for  tetragonal 
reduced  zone 

Co(l) 

Cr(2) 

Cr(8) 

Cr(4) 

Cr(6) 

Cr(5) 

*  IK-20  (72) 

±  0.53 

+  1.35 

+  0.97 

-0.41 

+  0.36 

-0.33 

TABLE  IV.  Spin  magnetic  moments  for  8-iayer  bcc  superlattices  Co(VIII)  and 
Co(III)Cr(V)  —  Co(l)Co(2)Co(3)Cr(4)  ...  Cr(8).  Atomic  sites  in  same  column  are  equiva¬ 
lent.  Note  two  distinct  solutions. 


Co(VHI)  < -  Co(III)Cr(V)  - > 


Mesh  for  tetragonal 
reduced  zone 

(ail  rites 
cobalt) 

Cod) 

Co(2) 

Co(8) 

Cr(3) 

CK7) 

Cr(4) 

Cr(6) 

Cr(5) 

IK-20  (72) 

IK-20  (72) 

+  1.72 
+  1.72 

+  1.S2 
+  1.53 

+  1.28 
+  1.38 

-0.29 
+  0.16 

+  0.22 
-0.22 

-0.11 
+  0.47 
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